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We extend our previous work on shot noise for entangled and spin polarized electrons in a beam- 
splitter geometry with spin-orbit (s-o) interaction in one of the incoming leads (lead 1). Besides 
accounting for both the Dresselhaus and the Rashba spin-orbit terms, we present general formulas for 
the shot noise of singlet and triplets states derived within the scattering approach. We determine 
the full scattering matrix of the system for the case of leads with two orbital channels coupled 
via weak s-o interactions inducing channel anticrossings. We show that this interband coupling 
coherently transfers electrons between the channels and gives rise to an additional modulation angle 
- dependent on both the Rashba and Dresselhaus interaction strengths - which allows for further 
independent coherent control of the electrons traversing the incoming leads. We derive explicit 
shot noise formulas for a variety of correlated pairs (e.g., Bell states) and lead spin polarizations. 
Interestingly, the singlet and each of the triplets defined along the quantization axis perpendicular 
to lead 1 (with the local s-o interaction) and in the plane of the beam splitter display distinctive 
shot noise for injection energies near the channel anticrossings; hence, one can tell apart all the 
triplets, in addition to the singlet, through noise measurements. We also find that spin-orbit induced 
backscattering within lead 1 reduces the visibility of the noise oscillations, due to the additional 
partition noise in this lead. Finally, we consider injection of two-particle wavepackets into leads 
with multiple discrete states and find that two-particle entanglement can still be observed via noise 
bunching and antibunching. 

PACS numbers: 71.70.Ej,72.70.+m,72.25.-b,73.23.-b,72.15.Gd 



I. INTRODUCTION 

Spin-related effects underlie promising possibilities 
in the emerging field of semiconductor spintron- 
ics and spin-based quantum computing^. Spin- 
entangled electron pairs in unconventional geometries, 
e.g., electron beam splitters^, offer a unique setting 
in which to investigate fundamental non-local elec- 
tron correlations in solids^. Several schemes for cre- 
ating and injecting entangled pairs in mcsoscopic sys- 
tems have recently been proposed involving quantum 
dots, superconductors, and interference in the electron 
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tion, coherent manipulation, and transfer of spin entan- 
glement ("flying qubits") in nanostructures are crucial 
ingredients for quantum-information processing and com- 
munication. Non-equilibrium noise, shot noise, is a useful 
probe for detecting entanglcmcn^i2 6 i. 

More recently, the Rashba spin-orbit interaction 
present in confined electron systems lacking structural 
inversion asymmetrjsSi has been proposed as a conve- 
nient means to spin rotate entangled pairs 2 ^. Interest- 
ingly, it was found that a local Rashba spin-orbit in- 
teraction acting upon a non-local portion of spatially- 
separated entangled electron pairs injected into a beam 
splitter gives rise to sizable modulation of the shot noise 
in the outgoing leads 2 ^. The use of the Rashba inter- 
action to controllably rotate the electron spin was first 
proposed by Datta and Das^ 9 -. Motivated by this ear- 
lier proposal and its potential impact on semiconductor 



spintronics, many researchers are actively investigating 
spin-orbit related physics in a variety of semiconductor 
nanostructureo 30 ! 31 ! 34 ! 35 ! 36 ! 37 ! 38 : 39 ! 40 ! 41 ! 42 ! 43 ! 44 ! 45 . 

Here we extend our previous investigation on the co- 
herent s-o control of entangled and spin-polarized elec- 
trons and their shot noise for transport in a beam-splitter 
configuration (Fig. 1) with local spin-orbit interactions, 
i.e., interactions acting within only a finite region of 
one of the two one-dimensional incoming leads 2 ^. We 
include both the Rashba 2 ^ and the Dresselhaus 4 ^ spin- 
orbit terms^ 7 -. Since the Rashba part of the s-o coupling 
is gate-tunable^ 8 ., one can controllably spin rotate the 
incoming correlated spinor pairs thus changing the de- 
gree of symmetry of the spin part of pair wave function. 
The stringent requirement of antisymmetry for fcrmions 
- the Pauli principle - intrinsically links the spin and 
the orbital (charge) degrees of freedom^. Thus the spin- 
orbit induced spin rotation affects the spatial charge dis- 
tribution of the pair which can be probed via current- 
fluctuation measurements: charge shot noise. 

We consider a beam-splitter with quasi-one dimen- 
sional incoming leads with one and two channels, (i) For 
singlc-moded leads and within the scattering approach 
we generalize our previous results 2 ^ by deriving general 
expressions for the shot noise of singlet and triplet pairs 
injected into the beam splitter. We present explicit for- 
mulas for the particular beam-splitter scattering matrix 
of the experiment in Ref . and a variety of incoming elec- 
tron pairs: singlet and entangled and unentangled triplet 
states defined along distinct quantization axis, (ii) The 
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FIG. 1: (a) Spin-entangled electrons injected into a beam- 
splitter setup with spin-orbit interactions, Rashba and Dres- 
selhaus, within a finite region L of lead 1. The strength a of 
the Rashba interaction can, in principle, be controlled via a 
top gate so as to be equal or unequal to the Dresselhaus cou- 
pling (3. For two orbital channels in lead 1 and a = P, no s-o 
induced band mixing occurs, right panel (b). For a 7^ j3 (or 
when either a = or fi — 0) the bands anti cross, left panel 
(b). Only a single spin rotation 8 BO = 2m ^/ a 2 + {3 2 L/h 2 is 
present for a = j3, while an additional "mixing" spin rotation 
6d modulates the electron transport in lead 1 for a ^ j3 and 
impinging energies near the crossing e ~ e c . This modulation 
appears in the current fluctuations (shot noise) measured in 
lead 3. In particular, each of the triplets - for a quantization 
axis along the y direction - exhibits a distinctive noise as a 
function of (6 so ,6d)- 



case with two channels is particularly interesting as the 
s-o terms give rise to inter-channel coupling which results 
in anticrossings of the bands. For incoming energies near 
these avoided crossings, we find similarly to Ref. an 
additional spin phase due to the coherent transfer of car- 
riers between the s-o coupled bands. Here, however, this 
modulation angle depends on both the Rashba and the 
Dresselhaus coupling strengths. Interestingly, for singlet 
and triplets defined along the y quantization axis (Fig. 
1) and injected into only one of the two channels of the 
incoming leads, we find that each of the triplet (besides 
the singlet) pairs displays distinctive noise modulations. 
This provides a way of distinguishing all of these triplet 
pairs via noise measurements. The intcrband coupling 
controlling the extra phase can, in principle, be varied 
via independent side gates which changes the width of the 
incoming channels^; this provides an additional mecha- 
nism for electric spin control. Moreover, for tuned s-o 
couplings (i.e., equal strengths) the Rashba and Dressel- 
haus terms partially cancel themselves out, thus giving 
rise to parabolic-band crossings for arbitrary strength of 
the s-o interaction^, Fig. 1 (b). This allows for the prop- 
agation of electron spins protected against non-magnetic 



scattering, i. e., robust entangled or unentangled spin 
pairs. 

We also consider spin-polarized inject ionS^i into the 
beamsplitter. Here we find that noise measurement can 
probe the spin-polarization of the Fermi liquid leads 
along distinct quantization directions. We also discuss 
the effects of backscattering in the incoming leads, due to, 
e.g., the potential discontinuities at the entrance and exit 
of the s-o active region in lead 1 (see Appendix iBl for an 
explicit evaluation of the transmission coefficient for elec- 
trons crossing a ID lead with s-o interaction). Backscat- 
tering reduces the visibility of the shot noise oscillations, 
because of the additional partition noise in the incom- 
ing leads. Finally, we investigate transport of injected 
two-particle wavepackets into leads with multiple discrete 
states but without s-o interaction. Similarly to our previ- 
ous results^ with ordinary two-particle pairs (i.e., "plane 
waves"), we find that two-particle entanglement can also 
be detected via noise measurements (bunching and anti- 
bunching) even with incoming wavepackets. 

This paper is organized as follows. In Sec. II we in- 
troduce the spin-orbit Hamiltonian in ID channels. We 
consider both the Rashba and the Dresselhaus s-o terms. 
We present exact and approximate solutions for wires 
with, respectively, equal and unequal (Rashba and Dres- 
selhaus) s-o coupling strengths. The full s-o transfer ma- 
trix for wires with one and two (coupled) channels is also 
derived. The boundary conditions for the two coupled 
channel case is discussed in detail. In Sec. Ill we present 
the basics of the scattering formalism for current and shot 
noise of spin-entangled electron pairs and spin-polarized 
electrons. Wc derive general formulas for the shot noise 
of singlet and triplet pairs injected into a beam splitter 
with an arbitrary scattering matrix (III-B) . The effect of 
backscattering is also discussed (III-C) for electron pairs 
in single-moded incoming leads. We present many spe- 
cific formulas for the noise of Bell pairs, electron pairs 
defined along distinct quantization axes for both singlc- 
and double-moded wires. Noise for spin-polarized injec- 
tion in discussed in (III-D). We also consider (III-E) the 
injection of entangled and unentangled wavepackets into 
leads with multiple energy levels. We summarize our re- 
sults and conclusions in Sec. IV. Many technical details 
of our calculation are discussed in the Appendices A-E. 



II. SPIN-ORBIT COUPLING IN ID 
CHANNELS: RASHBA AND DRESSELHAUS 

Quantum wires can be defined from two-dimensional 
electron gases by further constraining the electron motion 
to one spatial direction via, for instance, gate electrodes. 
When the underlying 2DEG has spin-orbit interactions 
of the Dresselhaus^ and Rashba 2 ^ types, due to bulk in- 
version asymmetry (BIA) and structural inversion asym- 
metry (SIA), respectively, the ID channel so formed will 
also present such interaction terms 02 . The Hamiltonian 
of a 2DEG with spin orbit interaction and an additional 
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gate-induced confining potential V(y) reads 

h 2 ( d 2 d 2 \ , , 

ia(a y d x — a x d y ) + ift(a y d y - a x d x ) , (1) 

where di = d/di, i = x,y and the third and fourth terms 
are the usual Rashba (strength a) and the linearized 
Drcsselhaus (strength ft) s-o terms, respectively. 

A. Exact solution: a = j3 case 

Similarly to the two-dimensional case treated in Ref. 
l49l the s-o wire problem here is exactly solvable for tuned 
couplings I a I = ft. Let us first consider the general 
case of a two-dimensional electron gas with an arbitrary 
scalar potential V(r) which can, e.g, describe static non- 
magnetic impurities, or further confinements creating a 
quantum wire or a quantum dot. At the symmetry points 
a = ±/3 the operator £ = (a x =F a v )/\2 provides an ad- 
ditional conserved quantity, and a general cigenstate of 
Tt and £ reads (for a — —ft) 



(2) 



in the o~ z basis. The function (p(r) fulfills the usual spin- 
independent Schrodinger equation 



and e is the energy eigenvalue of the wave function 
ip±(f) with £ = ±1. Now consider a quantum wire 
along the x-direction, i.e. V(r) = V(y). At a = —ft 
the wave functions are of the form (J2J with ip n (r) = 
4> n {y) exp(i(k ± y/2am/H 2 )x) / '\[L^, L x is a normalizing 
length, such that the full wave function reads 



V2 V ±eW4 



ik.r 



--<t>n{y)e 



±i y/2 OLiny / h 2 



(4) 



where <fi n {y) obeys the usual Schrodinger equation for 
the transverse variable y with quantized eigenvalues e n 
. The eigenstates Q are characterized by the subband 
index n and the wave number k, and the corresponding 
cigenencrgies are given by = s n + (h 2 /2m)(k ± 



\/2am/h 2 ) 2 — 2a 2 m/h 2 . Note that, similarly to the two- 
dimensional case discussed earlier—, the wire energy dis- 
persions here are also parabolic - for any strength of the 
\a\ = ft coupling, see Fig. 1(b). 



B. Approximate solutions: a//? case 

For unequal couplings we first solve the quantum wire 
problem in the absence of spin orbit coupling and then 
use this solution as a basis to write down the Hamilto- 
nian matrix with the s-o terms. Here we neglect any ad- 
ditional s-o terms arising from the further confinement 30 

V(y). 

1. Quantum wire eigenstates 
The solution to Eq. without the s-o terms is 



(p k ,n,*A x ,y) = -7f=<l>n(y)Wz), 



(5) 



where \<j z ) E {| |) z | is the electron spin state in the 
o~ z basis, with eigenvalues 



21.2 



£fc ' 



H 2 k 
2m 



(6) 



and n = a, 6... denoting the transverse modes with ener- 
gies e n (note that e n = e„ in the absence of s-o). The 
transverse confining eigenfunctions 4> n {y) obey the ID 
Schodinger equation 



h 2 d 2 Mv) 

2m dy 2 



+ V(y)<f> n (y) = e n </> n (y). 



(7) 



The confining potential in Eq. (JJJ is arbitrary. Later 
on we consider an explicit form (obtained for hard-wall 
confinement) so as to obtain simple estimates. 



2. Rashba-Dresselhaus wire 

We can derive a reduced Hamiltonian for our quan- 
tum wire with s-o by expanding the solution of 
Eq. Q in the basis of the wire without s-o, 

{<Pk,a,h<Pk,a,i> <PkM, <Pk,b,l}- Hcre we consider only two 
wire modes. Wc then find 



^-+e a (ia + ft)k (-ia + ft)d ab 

{-ia + ft)k ^- + e a (-ia-ft)d ab 

(ia-ft)d* ab ^L + e b (ia + ft)k 

(ia + ft)d* ab (-ia + ft)k ^ + e b 
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The matrix element 



with 



dab = -d* ba = (cj) a \d/dy\<f> b ) 



(9) 



in Eq. (JHJ) defines the s-o induced interband mixing be- 
tween the wire modes arising from the s-o terms pro- 
portional to p y in Eq. JQ). For hard- wall confinement 
d a b = 8/3w, where w is the wire width. It is convenient 
to rewrite the above matrix in the basis of the eigenstates 
corresponding to d a b = 0. For null interband coupling the 
Hamiltonian decouples into two sets of s-o bands 



2 J.2 



h 2 k 

2m 



+ e„ — sky' a 2 + (3 2 , 



(10) 



where n = a, b and s = ±, and eigenvectors 



¥>k,n, S { x ,y) 



1» 



4>n(y)\s), 



ti(-<) = ^ (IT> -- <u> ' ) ' (12) 



where 



f = Va 2 + P 2 / (ia + /3) 



(13) 



We define the transformed Hamiltonian matrix as H 
WHU, with 




(14) 



(11) We find 



H 







-2id 



a/3 





,* a 2 -0 2 

2id* , aP 

ab ^W 2 



2id a i 
d a t 



V" 2 +/3 2 
o?-0 2 



-d a l 



a 2 -p 2 

-2id, ■ 



ab ) — z w 

y/a 2 +P 2 





(15) 



The diagonalization of l|15(l is straightforward; the 
eigenenergies are 




where s, s' = ± . The corresponding cigenfunctions arc 
too lengthy to be shown here. Figure 1(b) shows the 
above energy dispersions for a ^ (3 and a = (3 for nonzero 
interband coupling d a b- In general, the energy dispersions 
present avoided crossings for a ^ (3. In contrast, the s-o 
tuned a — (3 case has eigenvalues which are quadratic in 
k with no avoided crossings. This k dependence is easily 
seen by setting a = (3 in Eq. I|16|) 



£s.s>{k) = - — + - (e b + e a ) + sV2ka + 
2m 2 

s'lV(tb-ea) 2 + 8a 2 \dab\ 2 . (17) 



In what follows we discuss in more detail the cases 
dab = and d a b ^ corresponding to the uncoupled and 
interband-coupled channels, respectively. We emphasize 
again that the interband coupling described by the ma- 
trix element d a b is purely induced by the s-o. As we 
will sec below, the uncoupled case gives rise to a single 
spin-rotation modulation. The interband coupled case, 
on the other hand, will have two independent modula- 
tion angles for injected electrons with energies near the 
band crossings. 
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Uncoupled ID channels (d a ij — 0): single spin 
rotation Or 



the absence of s-o induced channel coupling (uncoupled 
channels) by the 4x4 "transfer" matrix U" Q 



Here we have in mind a two-terminal geometry with 
the source and drain connected by a Rashba-Dresselhaus 
wire. For simplicity, we neglect the band offsets between 
the various interfaces. That is, we assume a unity trans- 
mission through the s-o region^. Finite offsets give rise 
to Fabry-Perot type oscillations which further modulate 
the transport properties^! of the system. The uncou- 
pled case (d a b = 0) considered here should be a good 
approximation also for finite d ab , provided that a\d a b\ 
be much smaller than the interband energy separation 
(a|d a h| -C £(, — e n ). The solution for d ab = is straight- 
forward (see Ref. |3l] for the case where only the Rashba 
coupling is active). From Eq. (|1 5|) . which is diagonal for 
d a b — 0, we immediately obtain the two sets of s-o bands 
[Eq. HlOJl ] which we rewrite as 



where 



11. 

2m 



(fc — sk so ) + e a ,b 



2 1,2 



K 2 k 



2m 



mV'a 2 +(3 2 /h 2 



s = ±, (18) 



(19) 



is the s-o wave vector. The corresponding eigenvectors 
are given in Eqs. I|llfl an d i|12|) . For d a b = the s-o 
bands cross at 



(-b - e a 
2V« 2 +/3 2 



(20) 



which is obtained by setting Ea(k c ) = e^(k c ) (see thin 
solid line in the inset of Fig. 2); a symmetric crossing 
also occurs for at k = — k c . 

As first pointed out by Datta and Dais 2 ^, injected elec- 
trons moving down the ID channel will spin precess due 
to the action of the s-o interaction. Here the spin ro- 
tation is due to the combined effects of the Rashba and 
Dresselhaus terms. In analogy to the case discussed by 
Datta and Das, here we find that a spin-up electron, say 
in channel a, crossing the length L of the s-o active region 
will emerge in the state 



where 



| T>, -> cos(0 so /2)| t), - sin(0 so /2)| i>*, (21) 



6 so = 2m^a 2 + P 2 L/h 2 (22) 



is the spin rotation angle about the y axis. Similarly, a 
spin down electron evolves into 

| l) z -> sin(0 so /2)| t)« + cos(0 so /2)| I),. (23) 

The same reasoning applies to impinging electrons in 
channel b. Hence, we can described the s-o region in 



U" 



u^ o 
ut 



where 



u a = u ft 

^ so w s 



cos(0 so /2) sin((9 so /2) 
- sin(0 so /2) cos(0 so /2) 



(24) 



(25) 



defines the single-channel transfer matrix for the uncou- 
pled channels a and b. Later on we introduce the scat- 
tering matrix approach to calculate current and noise in 
a beam-splitter geometry. The s-o rotation matrix above 
(and its generalization for two channels) will prove very 
convenient in accounting for s-o effects on the transport 
properties of the beam-splitter within the scattering ap- 
proach. Note that only the Rashba coupling constant 
appearing in the rotation angle 9 so can be varied exter- 
nally via a gate electrode, while the Dresselhaus coupling 
P is a material property. As a final point, we note that 
the above s-o rotated states satisfy the proper boundary 
conditions for the wavefunction at x = and x = L. 
This is discussed in some detail in appendix A for both 
the one- and two-channel cases. 



D. Coupled ID channels (d a b 7^ 0): additional spin 
rotation f5 d for a/(3 

For nonzero s-o induced interband coupling d a bi the 
subbands anticross for distinct coupling strengths a 7^ 
p. Similarly to the one-channel case, here we also have 
to find out how incoming spin up (or down) electrons 
emerge after traversing the s-o active region of length L. 
Here we have in mind incoming electrons with energies 
near the d ab — crossing of the bands at k c , i.e., e ~ 
e~(k c ) = £i~(fe c ). This is the relevant energy range where 
s-o induced interband crossing should play a role (unless 
a = P). In what follows we present a simple analysis of 
this injection problem by using a perturbative approach 
("near free electron model"—) to describe the s-o states 
near the crossings. 

For injection energies near the d ab = crossings, we 
can approximate the Hamiltonian in Eq. I|15f) by 





d 
V 





< 

a 2 -j3 2 









d a \ 







\ 




4 J 



(26) 



i.e., we drop all the off-diagonal matrix elements except 
those directly coupling the states near the crossing. From 
the form of H app it is obvious that the crossing states 
[middle block of Eq. Q26[l] will split due to the d a b cou- 
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pling. The new eigenvalues are 

^27,2 



e±(k) 



where 



^k 2 1, , a 2 -/3 2 
o — + o ( £ b + e a ) ± \d ab \ = 

2"l 2 yV + ^2 



(e b - e a ) - 2 v /« 2 +/3 2 fc 



4 Kb 



(27) 
(28) 



can be viewed as an expansion parameter near k c [Eq. 
I|2()|l]. Expanding e±(k) near k c (we should keep only 
the lowest order in a; since the third term of Eq. (|27fl is 
already proportional to d ab ), we find to zeroth order in x 



e±(k) 



P 2 



2m 2 ^ a 2 + £2 



The corresponding eigenvectors are 
1 



V2 



(|->.±l+>&), 



(29) 



(30) 



where |-) Q -> ip k , a ,-(x,y) and |+) 6 -» (^^(x, y) are 
the eigenstates in Eq. (fTTf) - The new eigenstates \ip±) arc 
zcroth-order linear combinations of the crossing states 
(remember that the energies are linear in ad a b). More 
explicitly, 



1-0=1 



-(' 

2U 



(31) 



In a "four-vector notation" we can write 




(32) 



For small s-o induced intcrband coupling we look for sym- 
metric solutions around k c [Eq. (|20|) ]: fc cl = k c — A/2 
and k C 2 — k c + A/2. Equation l|33(l then gives 



2m\d a b 



a 2 



k c fr 2 v'a 2 + /3 2 ' 



(34) 



Having determined the wave vectors fc c i and fc C 2, we can 
now solve the injection problem. The idea is to expand 
the incoming electron state, say spin up in channel a, in 
terms of the eigenstates of the s-o region. The expansion 
has to satisfy the boundary conditions (continuity of the 
wavefunction and flux conservation) at both the entrance 
and the exit of the s-o region. 




k cl k c2 k 2 



FIG. 2: Schematic of the quantum wire energy dispersions 
£ s,s'{k) (Eq. 1161 for a/(3. The blowup shows the band anti- 
crossing for d 7^ in more detail. The crossing thin solid lines 
represent the uncoupled case, d a ^ = 0. The curves with circles 
are obtained from Eq. 1291 (e±(k)) and are good approxima- 
tion for the actual dispersions near crossing point The 
wave- vectors k c i, k C 2, and ki, used to expand an incoming 
plane wave within the s-o region [Eq. 1351 1. are also shown 
in the inset. 



Boundary conditions 



As Fig. 2 clearly shows, s±(k) [Eq. I|29|l ] approximate 
well the exact energy dispersions e SyS t (k) [Eq. of the 
problem near k c . By using Eq. (|29|) we can analytically 
determine the wave vectors k c \ and k C 2 relevant for the 
spin injection problem. This is easily done by imposing 
Sf = £+(k c i) = £-(k C 2) which yields 



&k 2 c2 
2m 



2m 



2\d ab 



P 2 



(33) 



Here we show that spin injection with energies near 
the band anticrossing is possible in our system, provided 
that the s-o interband coupling be small compared to the 
Fermi energy. Details are given in Appendix A. 

a. Continuity of the wave function. A spin-up elec- 
tron in channel a entering the s-o region at x = with 
an energy ef ~ £+{k c i) = £-(k C 2) has to satisfy 






Vo 




Ak c \x 



l 

4 



1 



-1 | 




//,_/./• 



(35) 



The above condition is clearly fulfilled; a similar condi- b. Continuity of the current flow. The continuity of 

tion holds at x = L (see Appendix A). the (non-diagonal) velocity operator— acting on the wave 
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function at x — which assures current conservation yields 



m 






ikx I 

e U-»o- 



/ ^- (fc c — A/2 + k so ) \ 
M (fc c _ A/2 + fc so ) 
A (fc c - A/2 - fc so ) 

V -e£ (fc c - a/2 - fc so ) / 




. /7,-_>.f 



/ A (fc c + A/2 + fc so ) \ 
(fc c + A/2 + k so ) 
-l(fc c + A/2-fc so ) 
V ( fcc + A/2 - k so ) J 



(36) 



which simplifies to 





m 




( 

1 
















~ 2 


I 





) 


V 



£(fc c +fc 2 ) 

^ (fcc - fc 2 - 2fc s 
-AA/2 
PA/2 





( 


(kc+fe>) 
2k F 


\ 


hkp 













A 




m 




4k F 




J 


{ 




) 



(37) 



where we have used fc 2 — fc c = 2fc so [Eq. I|18|) ]. From Eq. 
(|37|l we see that the matching of the derivative is fulfilled 
provided that A <C Akp. As we show later on, this is 
the case for realistic parameters. The velocity operator 
matching at x = L holds similarly (see Appendix A). 



2. General spin-rotated state at x = L 

After traversing the s-o region a, say, spin-up electron 
in channel a is described by the state 




Straightforward manipulations lead to 



\U. — _ p i(ka+k so )L 



/ cos (0 d /2) e - t9 °°/ 2 + e l9 °°/ 2 \ 
€ [cos (0 d /2) e -««/ 2 - e""/ 2 ] 
-isin(0 d /2)e-^°/ 2 
j£sin (0 d /2) e -*"/ a 

(39) 



/ 



where we have introduced the additional modulation an- 
gle 

9 d = AL= (\d ab \/k c ) 6 so (a 2 ~ (3 2 )/(a 2 + (3 2 ) (40) 

due to s-o induced interband mixing. We show in the 
Appendix A that for x ^ L, the state 



cos {6 d /2)e~ l9 »°/ 2 H 
e[cos(V2)e^°/ 2 



6so/2 

»e 3 o/2i 



i e *(*.+*")*0 o ( v ) 



-ism((V2)e- ies °/ 2 
^sin(6» d /2)e- ie -/ 2 



- fc "> a &(tf), (41) 



satisfies the proper boundary condition for the velocity fulfilling the continuity of the wave function at this in- 
operator [note that setting x = L in (|4*TJl gives i|3*9"|) , thus 
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terface]. Equation (|41|l shows that upon traversing the 
s-o active region of length L, a spin-up electron in the in- 
coming channel a, acquires a spin-down component in the 
same channel and, more importantly, coherently transfers 
into channel b. This coherent transfer from channel a to 
channel b is solely due to the s-o induced interband cou- 
pling near fc c , described by the mixing angle 9 d . Hence, a 
weak s-o induced interchannel mixing - rather than be- 
ing detrimental to transport - offers a unique possibility 
for further spin modulating the electron flow. 



feet of the s-o interaction on electrons impinging near 
the band crossing at k c . This transfer matrix is readily 
contructed in terms of the column vectors similar to the 
one in Eq. I|39|l . which describes how a spin up electron 
in channel a evolves upon crossing the s-o region. We 
obtain 



3. s-o transfer matrix: coupled channels 



Similarly to the uncoupled-channel case, here we can 
also define a s-o transfer matrix Ug describing the ef- 



U' 



2 6 



ik c L 



( cos (0 d /2) + e ie - 
£ [cos (6 d /2) - e u 
-i sin (0 d /2) 
\ i£sin(0 d /2) 



C [cos(0 d /2) " 
cos(<9 d /2) + e l9 °° 
-if sin {9 d /2) 
ism{9 d /2) 



-isin(<V2) 
-ifsin(0 d /2) 
cos(0 d /2) + e~ ie > 
-i [cos(0 d / 2 ) - e 



-ikd s 



iC sin(0 d /2) 
* sin (p d /2) 
-C* [cos(^/2)-e- 
cos (6 d /2) + e- ie * 



(42) 



where the modulation angles 9 so and 9 d are given in Eqs. 
()22J) and (|4()(l . respectively. We should keep in mind 
that Eq. I|42|) describes electrons traversing the s-o re- 
gion with energies near the crossing energy. As we dis- 
cuss later on, the s-o transfer matrix above is also useful 
for spin-rotating entangled and/or unentangled electron 
pairs injected into a four-terminal geometry (beam split- 
ter). The idea is that U cc operates on the member of 
the pair traversing the s-o region. Note that the transfer 
matrix in Eq. fll^ reduces to that of the uncoupled case, 
Eq. for a = /3 (9 d = 0). Next we estimate the 

magnitude of the spin rotations we have described here. 



E. Single spin rotation 3O for coupled channels 
(d ab / 0) with a = f3 

Here the calculation is simpler since the bands do not 
anticross even for non-zero d a b as we discussed earlier. 
The crossing wave vector k c for a = (3 is determined from 
Eq. (|r7Jl . For instance, the k > crossing is obtained by 
setting e +! _(fc c ) = e_ i+ (fc c ) which gives 

k c = -4- V(tb-e a ) 2 +8a 2 \d ab \ 2 . (43) 
2v2a 

For dab = and a = (3 the above wave vector reduces to 
k c defined in Eq. ||2UJ|. 

By expanding the incoming electron states into the ex- 
act eigenstates derived in Sec. IIA we can obtain the 
modulation angle 9 so = 2\/2mL/h 2 . Note that 9 d = 



for a. = j3. Interestingly, the matching of the bound- 
ary conditions here and the general state at x — L can 
be straightforwardly obtained from the a / |3 case by 
setting A — (or cquivalently, 9 d = 0). However, it is 
important to note that the crossing wave- vector is now k c 
(not k c ) and that k^ — k c = 2k so , where k so is calculated 
for a = (3. Note also that only one modulation angle 9 so 
is present for the tuned-coupling case a = (3. Hence, this 
case is similar to the uncoupled channel problem treated 
by Datta and Das^, even though here d a b ^ 0. The 
identical coupling strengths makes the problem similar 
to that of the uncoupled channels; however the rotation 
angle is now renormalized. 



F. Estimates for the modulating angles 6 so and 6a 

Simple estimates for the spin-rotation angle 9 so and the 
mixing angle 9 d can be obtained by assuming a hard-wall 
transverse confinement of width w. Using the well-known 
analytical results for the wire problem, we find d a b = 
8/3w for the intcrband mixing and eb = 3ir 2 h 2 /2mw 2 
(assuming e a = 0). The quantity e so = h 2 k 2 /2m = 
m (a 2 + /3 2 ) /2h 2 sets an energy scale in our problem. 
For the sake of concreteness, let us choose Cb = 16e so 
which leads to y 'a 2 + (3 2 = {^/S/2n / '4) ft 2 / 'mw = 2.44 x 
10~ 2 eVnm (and e so ~ 0.2 meV) for m = 0.05mo (see 
Ref. [48^ and w = 60 nm. The energy at the band 
crossing points is then e~(k c ) = 24e so ~ 4. 8 meV; 
note that for Fermi energies close to this value, s-o in- 
duced channel- mixing effects are important. From Eq. 
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(I20|l we find k c = 8e so / \/ a 2 + (3 2 . Assuming an ac- 
tive s-o region of length L = 69 nm we can estimate 
the spin-rotation angles; we find 9 so = tt. To obtain 
6d = (dab/kc) 9 so (a 2 — /3 2 )/(a 2 + j3 2 ) we need an esti- 
mate for /?. To estimate the Dresselhaus coefficient in 
a quantum well geometry we use /? = ^(k 2 ), where (k 2 ) 
denotes the expectation value of the wave vector compo- 
nent along the growth direction. For the lowest infinite- 
well eigenstate we find {k 2 ) = (ir/w) 2 . The coefficient 

7 is typically w 25 cVA 356 i 57 i 58 which yields (3 w 1CT 5 
eVnm. Hence, for such III-V materials we can neglect 
(3 and use 9d = {d a b/k c )9 so , which gives 9d = tt/2 since 
dab/K = 2/(3k so w) ~ 0.5. 

In order to obtain comparable Rashba and Dresselhaus 
coupling strengths, we could use a setup with wider wires 
and materials with a larger effective mass^S. In addition, 
we could consider an inhomogencous beam-splitter with 
a different material with larger Dresselhaus coupling in 
one of the incoming arms. Note that the possibility of 
tuned couplings a = [3 is very attractive since in this case 
the spin of the electron propagating in the s-o coupled 
channels is insensitive to non-magnetic impurity scatter- 
ing HII Ap . i.e., the spinor is k independent for a = j3. 

We stress that the modulation angles 9 so and 9d can, 
in principle, be tuned independently via a proper gat- 
ing structure. This could involve, for instance, both side 
(top) and back gates to induce changes in the channel 
width w (confining potential) and the Rashba constant. 
The above conservative estimates suggests that the spin 
rotations we are considering here are sizable. Finally, we 
note that for the above parameters A/Akp ~ 0.05 <C 1, 
which justifies the approximation made in the velocity 
operator matching [Eq. (|37|) ]. 



III. TRANPORT PROPERTIES: CURRENT 
AND NOISE 

In what follows we calculate the current and its dy- 
namic fluctuations (shot noise) for electrons traversing 
a beam-splitter. We use the scattering approach of 
Landauer and Buttiker— . We consider injection of (i) 
electron pairs (singlet and triplets) from an "entangler" 
tunnel-coupled to the incoming leads of the beam-splitter 
and (ii) spin-polarized electrons from Fermi-liquid leads 
which arc assumed to be thermal reservoirs each held at 
a given chemical potential. For a calculation of shot noise 
for entangled electrons in a beam-splitter where a Berry 
phase provides an additional modulation, see Ref. Iriol 



A. Scattering approach: basics 

Here we briefly outline the scattering-matrix formula- 
tion for current and noised. 



1 . Current 

Within the Landauer-Buttiker approach, the transport 
properties of a mesoscopic system are expressed in terms 
of the scattering matrix s w connecting the many incom- 
ing and outgoing attached leads. The current operator 
in lead 7 is 

t r( t ) = TZ E E A % (r, e, e')e^-*'WaUeW«> (O, 

a.paa' es 

(44) 

with 

(45) 

where a =|,| is the relevant spin component along 
a proper quantization direction ( u x, y, or 2"). Wc 
have introduced the creation (annihilation) fermionic op- 
erator a^ acr (e) [a aa (e)] for an electron with energy e 
in lead a, which satisfy the anticommutation relation 
{al, a (e),a a 'a->(e')} = 6 aa >5 aa ,6 £e > We have considered 
beam-splitter leads with discrete longitudinal energy lev- 
els e,e'. This yields the factor v = {L x /2-kK) \Jra/2Ep 
in (|44|1 . which actually is the ID density of states for 
only forward propagating states (positive momenta). In 
the standard expression for the current with continuous 
energies^, this factor cancels with the density of states 
appearing when transforming discrete states into contin- 
uous ones. For a study of noise in a beam-splitter with 
continuous energies, see Refs. l6ll andl62l We discuss in 
more details the transition from the discrete case to the 
continuous one in Sec. IIII El 

2. Shot noise 

At a time t, the current fluctuation about its average in 
lead 7 is SI 7 (t) = I 7 (t) — (Iy). In a multiple-lead configu- 
ration, the shot noise between leads 7 and fi is defined as 
the Fourier transform of the symmetrized current-current 
autocorrelation function 

■M<") = 5 J(6^t)5iptf) + 5iptf)5i^t))e™db. (46) 

The angle brackets in Eq. I)46|l stand for either (i) a quan- 
tum mechanical expectation value between two-particle 
states or (ii) a standard ensemble average (thermal reser- 
voirs). Note that the non-equilibrium current noise de- 
fined above arises physically from the discrete nature of 
the charge flow in the system. This is strictly true only 
at zero temperatures; at finite temperatures Eq. I|46|) 
contains also thermal noise. 



3. Beam-splitter scattering matrix 

To calculate the noise from Eq. (|46|l we need to specify 
the beam splitter scattering matrix. For a symmetric 
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beam splitter without s-o interaction and single-mode 
channels, we have the scattering matris 3 ^ 



Shot noise for electron pairs: singlet and triplet 
states 



V 









Sl3 


Sl4 \ 




f° 





r 










S23 


«24 










t 


;) 


S31 


S32 





1 


- 


r 


t 







S41 


S42 





J 






r 





0/ 



(47) 



that is, the beam splitter transmits electrons between 
leads 1 and 3 and leads 2 and 4 with amplitude r and be- 
tween leads 2 and 3 and leads 2 and 4 with amplitude i. 
Note that backscattering is neglected in s; see Sec. III-C 
and Appendix IC II for a beamsplitter including backscat- 
tering effects. If the incoming or outgoing leads have 
more than one mode (i.e., many quantized channels) we 
can, in a first approximation, assume that the beamsplit- 
ter does not mix the orbital channels so that Eq. I|47|) 
holds true for each of the modes separately. 

Interestingly, in the presence of s-o interaction in lead 
1, Fig. 1, we can define an extended beam-splitter scat- 
tering matrix to incorporate the spin rotation described 
by the s-o transfer matrix U so . Since an electron in lead 
1 undergoes a spin rotation described by U so , we can 
redefine the matrix elements S13 = S31 and S14 = S41 as 
4x4 matrices 



and 



s 31 = s 13U so , 



smU so , 



(48) 



(49) 



to incorporate the effects of the s-o interaction. Note that 
U so is given by Eqs. (|24|) and (|42|) for the uncoupled- and 
the coupled two-channel cases, respectively. The other 
elements in s remain unaltered except that they are now 
4x4 matrices, e.g., S24 = S42 = S14I, where 1 denotes 
the 4x4 unit matrix. Note that the new beam splitter 
scattering matrix s so incorporating the s-o effects in lead 
1 as defined above is a 16 x 16 object as opposed to the 
4x4 matrix in Eq. P7|> . 



We assume that an 

cn ^ an „^ r 5.6. 7.8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. 24. 25 j g 

placed just before leads 1 and 2, Fig. l(a)S. Below we 
calculate the noise for the states 



\S) 
\T F . ) 



= -^K T (ei)a2j.( £ 2) T a^(£i)4 T (£ 2 ) 
= aL( £l )aL( £2 )|0), <7=U, (50) 



where |0) denotes the ground state (filled) Fermi sea 
of the leads and i = x,y, z any particular quantization 
axis. The states |5*) and \T ei ) arc entangled singlet and 
triplet, respectively, while |T UT .) and \T Ul .) are unentan- 
gled triplets. Here we consider zero temperatures, zero 
applied voltages and zero frequencies. In this limit the 
Fermi sea is completely inert (noiseless) and the noise in 
the system is solely due to the injected pairs above the 
Fermi surface^. To determine the shot noise we essen- 
tially evaluate matrix elements of the general form 

(0\a^ (e^K,^ (e u )al a (e)ap^ (e') x 
4',."(£")«/3',^(£"')4^( £ 7)4,^(^)|0>, (51) 

appearing in the noise definition, Eq. (|46|l . This is most 
systematically done via Wick's theorem since the object 
in (|48|) resembles a four-particle Green function (see Ap- 
pendix D for details). 



1. General noise formulas: single- channel case 

For the injected singlet and the triplets in Eq. i|50|) 
we find the following expressions for the zero- frequency 
noise between leads 7 and \i 
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- T„ 



e 

2hv 



and 

„2 



a 7M — ■ 



e 

hv 



a=l, 2,0=1.. 4,cr',<7 a/ / 3=l,2,cr 

a^/3=l,2;cr 

+ \ A a%{T^ a ,e a )A a f' y {ii]e l 3,e l 3) - - ^ A£*( 7 ; e a , e a )A^(n; £(3 ) 



q^/3=1 i 2;it 



a,/3=l,2;<7 



o ^Q,a(7; £ o,£a)^ a ,a ° (Mi £ a> e a) 



a=l,2:<7 



(52) 
(53) 



a=l, 2,0=1.. 4,cr' 



a,/3=l,2 



Note that Eqs. I|52[) and (|54|l do not depend on the par- 
ticular form i|47|) of the beam-splitter scattering matrix 
and the quantization axis chosen. In what follows we 
present explicit formulas for the noise derived from Eqs. 
(|52[l and 154|) . We also determine the noise for entangled 
and unentangled states defined along distinct quantiza- 
tion axes (i = x, y and z) and for the Bell states. Later 
on we present similar results for the two-channel case as 
well. 



a. Quantization axis along z. For the specific forms 



and 



'13 



Si A 



rcos(0 so /2) rsin(0 so /2)\ 
-rsin(0 so /2) rcos((9 so /2) J 



tcos{6 so /2) isin(0 so /2) 
-fsin(0 flO /2) tcos(e so /2) 



(55) 



(56) 



2. Specific formulas: uncoupled- channel case 

In the absence of s-o induced interband coupling 
(uncouplcd-channcl case), the channels a and b are in- 
dependent within lead 1. That is, if electrons are in- 
jected only in the channel a of lead 1, they will remain 
in that channel while propagating through the length L 
of the s-o region in that lead. If fact, the channel in- 
dex remains unaltered as the electrons traverse the beam 
splitter since we assume the the beamsplitter does not 
mix the channels. We present below results for electron 
pairs injected only in channel a of the incoming leads. 
The case with two pairs injected into channels a and b 
is straighforward (factor of two) since no s-o interband 
mixing is considered here. However, as we discuss later 
on, in the coupled-channel case injection into just one of 
the channels is significantly different from injection into 
both channels (not just a factor of two as here). Below 
we detail the calculation of the noise from Eqs. I|52|) and 



l)54p. To calculate the noise in lead 3, i.e., S 



S/T e . 



33 



and 



1 , we first have to determine the relevant elements 
a g (3;s a ,£ a ) [Eq. (|45J) ] appearing in these quantities. 



°33 
A a,a 



Since our scattering matrix is assumed to be independent 
of the energy, so is (3; e a , e a ) = (3). 



the only non-zero » (3)'s are 

4:1(3) 



., 2 (3) 
. 2 (3) 



A T < T 



-S 31;iT S32;lT - s 3l;TT s 32;TT 
-r*tcos(0 so /2) = A[] l 2 (3) 

~ S 31;1T S 32;11 _ S 31;TT S32 ;Ti 

r*t sin(0 so /2) = -Aj;|(3) 

- S 32;1T S 32;1T ~ S 32;TT S 32;TT 

i.a(3)- 



-\t\ 2 = A l ~ l < 



(57) 
(58) 
(59) 



Note that in the above we have chosen the a index in s so 
to be that of the z component of the spin: a — > a z =f 
,|, i.e., we have set the quantization axis to be z, Fig. 
1. Hence, the entangled and non-entangled triplet states 

here refer to this basis: S 33 ' and S 33 "' z ; the noise for the 
singlet state S 33 is the same for all quantization axes. 

Plugging in the above A^Z (3)'s into Eq. (|52|) we find 
2e 2 RT, 



S. 



33 



(*„,) 



hv 



-[l±cos(0 so /2) cos(0 so /2)<5 ei , 6 



sin(0 so /2) sin(0 so /2)<S eije 



(60) 



where we have defined the transmission and reflection 
probabilities T = \t\ 2 and R = \r\ 2 , respectively. Since 
R + T = 1 , further simplifications lead to 

, 2e 2 



S- 



33V^so 



hv 



T(l-T) [1 + cos {6 so )5 eu 



(61) 
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S. 



33 



^ 2 e 2 , 



r)(i-<y eilBa ). 



(62) 



Similarly, we find 



^33 '" 1 



= ^33 (Geo) 



2f_ 
hv 



T(l-T)[l-cos 2 (9 so /2)S £u£2 l63) 



for the unentangled triplets with spin polarization along 
z. 

The above formulas have been derived for the case 
where the injected electrons in arm 1 and 2 have the 
infinitely-sharp energies e\ and £2. In Sec. IIII El wc 
consider the case where the injected electrons are de- 
scribed by Lorcntzian wave-packets of width 7 centered 
on e\ and £2. We then find, in the continuous limit, 
that £ £1 , E2 in Eqs. ((61163(1 is replaced by the function 
H(A) =V/(A 2 /4 + 7 2) w here A = £1 - £ 2 . This func- 
tion, which interpolates between H = 1 when A -C 7 and 
H = when A 3> 7, corresponds to the one appearing 
in Eq. (28) of Ref. |M 

b. Quantization axis along y. To obtain the corre- 
sponding noise espressions for injected electron pairs with 
spin polarization defined along the y we have to first 
rewrite s so in the basis of a y : | j) s = (| T)z + i\ l)z) /V%, 
I i)y = (I T)» _ A i)z) /V%- For instance, sfg becomes 



'13,y 



rexp(z6» so /2) 

rexp(-*0 8O /2) 



(64) 



Now the only non-zero elements are A['2(3) 
— r*texp(— i9 so /2) = 



4'i(3) and Ai;{(3) 



T.T/ 



— rt* cxp(—i6 so /2) 



4i(3) 



Substituting these 



terms into the general Eqs. ((52(1 and 1(54(1 . we find 

2e 2 



SH T ' V (Oso) = j^T(l - T) [1 ± cos (0 SO ) * ei , e „] , (65) 



and 



^33 T " (@so) 



} x>> 



(66) 

c. i?e/Z states. For completeness we also calculate 
the noise for injected Bell pairsS^( maximally entangled 
states) described by 



The Bell states above are defined with respect to an ar- 
bitrary quantization axis. The noise expressions for the 
first two Bell states are the same as those of the sin- 
glet and the entangled triplet derived above. The noise 
for the states | \&2) and ^3) can be easily determined in 
terms of the results for the unentangled triplet states for 
particular quantization axes. Let us consider the quan- 
tization axis along z, for concreteness. The procedure is 
straightforward: 

i) After traversing the s-o region, the injected unentan- 
gled triplet state |T„ T J becomes \T U ^ Z ) L = U° |T„ T J 
and can be decomposed as 

|T„ T z ) L = [cos(0 so /2)| T)ia,*-sin(0 so /2)| l)i ,,]®| T>2a.z, 

(71) 

or 



cos(0 so /2)| TT)c 



-sin(0 so /2)| |T)aa,,, (72) 



where we have used the shorthand notation | tT)oo,« = 
I T)ia,z <8> I ])2a,z = \T U] to denote the tensor product 
state with one electron injected into the channel a of lead 
1 and another in the channel a of lead 2. Similarly, 



IT,, 



sin(0 so /2)| Tl) 



■cos(0 so /2)| U> c 



(73) 



ii) Now, we rewrite the rotated Bell state \^2) L explicitly 



1*5 



COS 



~T" I J, J,) em, 2: ) 



sm 



,/2)-^(|U>«a,,-|it>aa,»).(74) 



iii) The noise in lead 3 corresponding to \^2) L is deter- 
mined from the expectation value of Eq. ((46(1 in this 
state. Since the beam-splitter scattering matrix does not 
include spin-flip processes, this expectation value can be 
expressed in terms of the shot noise for the unrotated 
singlet (| U) aa ,z - I 4T)oo,z)/2 and unentangled triplets 

TT)cm,Z7 I -14) aa,z 
$33 (@so) 



l - cos 2 {e so /2)[sl;^ a 



s. 



33;aa J 



/2)^ 3 

; a a ■ 



where 



and 



OS 

J 33;aa 



hv 



T(l-T)(l + 5 e 



(75) 



(76) 



1*2) 
1*3) 



\S) 

\T ei 
1 

71 
1 

71 



(\T uhl ) + \T uU )) 
{\T u i,i) - \Tui,i)) 



(67) 
(68) 

(69) 
(70) 



S. 



33 :aa 



J 33 ;aa 



2e 2 

= — T(l-T)(l 
hv 



(77) 



Note that S§ s . aa and S 3 ^." a ^ denote the shot noise for the 
bare singlet and triplet states, respectively^. Equations 
((76(1 and 1(77(1 represent the noise expectation value for 
the bare singlet and triplet states^ and follows from Eqs. 
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and with 9 so = 0. Substituting the Eqs. CEJ 
and |77J> into f75)l. we find 

Similarly, 

^3* 3 (*«0 - — T) (1 — S EUE2 ). (79) 

It is interesting to note that = S 3 ^ v and Sf^ = 
S^ v = Sl; iv - This follows from \T Ey ) = -*-(| U}, - 
| J.1%) being equal to |* 2 ) and |* 3 ) = |T„ T J for the 
quantization axis along z. 

d. Local Zeeman and spin- orbit induced rotations. 
It is interesting to note that the shot noise modulation 
induced by the local s-o interaction in the uncoupled- 
channcl case discussed above is formally identical to that 
due to a local Zceman-interaction in lead 1 [e.g.. due 
to a local magnetic field or a distinct g factors in lead 
1]. For a local Zeeman interaction, 9 so should be equal 
to gfiBBL/v, i.e., the phase acquired by the electron 
upon traversing, with velocity v, the length L of the 
Zeeman-active region. Note, however, that this formal 
correspondence between local Zeeman and s-o interac- 
tions holds only for the uncouplcd-channel case (or in 
a strictly single-channel lead) and when magnctic-field- 
induced orbital effects are neglected. 

e. Shot noise for electron pairs: physical picture. 
Before moving over to the more elaborate case with cou- 
pled channels, we provide here a simple picture for the 
shot noise results we obtained for the singlet and triplets. 
Let us first consider the case with no spin orbit (9 so = 0) 
and one orbital channel. Our formulas yield bunching 
for singlet and anti-bunching for the triplets, as previ- 
ously found in Ref. 0. This bunching and antibunching 
behavior for the singlet and triplets pairs can be read- 
ily understood if we write out the pair wave-functions 
in the outgoing leads (3,4) similarly to Ref. |6^ The 
injected singlet and entangled triplets states \S/T e ,) = 
(a| )T 4,x T a{ tl alj)\Q) /y/2 and \T Ua i ) = a\ a al a \0) after 
trasversing the beam-splitter evolve into 

\S) y/2rtbl A bljO) + V2rtb\ A bljO) + 



I^J-^^-^^Xi + ^Xt)! )' ( 81 ) 

and 

irO-^-r'RX*. ( 82 ) 

where the b\ a 's (i = 3, 4) represent the creation operators 
in the outgoing leads, directly related to the aj >cr 's via the 



beam-splitter scattering matrix s (note that here the par- 
ticular quantization axis is irrelevant; this is not true in 
the presence of the spin-orbit interaction as described be- 
low). Interestingly, the incoming electrons in the triplet 
states in leads 1 and 2 have zero probability of emerging 
in the same outgoing lead (note that | r 2 — i 2 1 = 1 due to 
s^s = 1). This is the (full) antibunching we mentioned 
above and it results in zero shot noise for the injected 
triplet pairs. Physically, this vanishing of the shot noise 
means that there is no randomness in the electron flow in 
the outgoing leads: each electron of the incoming triplet 
pairs goes - with unity probability - to distinct outgo- 
ing leads. For the singlet pairs, on the other hand, the 
probability for the two electrons to emerge in the same 
outgoing arm is not zero as for the triplets. This proba- 
bility is larger by a factor of two than the classical value 
|r| 2 |t| 2 . Hence the incoming electrons in a singlet pair 
tend to bunch, i.e., to go to the same outgoing leads. 
Full bunching occurs only for the particular case of a 
50:50 beam splitter. In this case we have r = ±it and 
opposite spins have a 1/2 probability to emerge in one of 
the leads 3 or 4. This randomness in the electron flow in 
the outgoing leads increases the shot noise as compared 
to both the triplet and the classical case. 



Note that the bunching and anti-bunching behaviors 
described above follow from the stringent requirement for 
anti-symmetry (Pauli's principle) of the total wave func- 
tion of the electron pair. In an uncntanglcd triplet, for 
instance, the spin part of the wave function is symmetric 
thus forcing its spatial component to be anti-symmetric 
which results in a strong correlation between the elec- 
trons in the triplet pair, i.e., they avoid each other by al- 
ways going into distinct outgoing leads. The same is true 
for the entangled triplet pair. This correlation reduces 
the shot noise, in agreement with the classical notion of 
shot noise suppression which occurs when the "discrete- 
ness" of the electron flow is reduced. For the singlet, on 
the other hand, the spin part of the pair wave function 
is anti-symmetric which makes the spatial part symmet- 
ric. The electrons in a singlet pair can overlap freely in 
space and hence are less correlated than the triplet case 
thus giving rise to non-zero shot noise. Interestingly, the 
shot noise for the singlet is even larger than the classi- 
cal analog (i.e., a distinguishable pair injected into leads 
1 and 2). This also agrees with the classical notion of 
shot noise as due to the discreteness of the charge flow: 
because the spatial part of the singlet pair is maximally 
"uncorrelated" ("negative correlation") the discreteness 
of the electron charge is larger than the classical case (the 
singlet electrons can lie completely on top of each other 
in real space) thus yielding a larger shot noise. 



In the presence of spin-orbit coupling, the above ideas 
can be straightforwardly generalized. In this case the sin- 
glet and triplets injected into single-channel leads evolve 
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into 



evolves into [Eq. l(39"j) ] 
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b lAi 



blA, 



upon traversing the s-o region in lead 1 (we dropped 
an overall phase e i ( fe <=+ fc =°) L . Similarly to the uncouplcd- 

|0) + channel case, we note that only the portion of the electron 
pair going through lead 1 is subject to the s-o effect. Note 

10(584) that the s-o induced interband coupling in lead 1 makes 
the injected states initially in channel a leak into channel 
b. The expectation value of the noise in the state |T„ T z ) L 
is 



\T UUZ ) cos(# so /2)(i 2 -r 2 )^ T ^ T |0) + 

sin(0 so /2)rt (4, T 4,l + b lA,l) 1°) + 
sin(0 so /2)t 2 &+ !T &ijO) - 
sin(^ o /2)r 2 6^6t T |0). (85) 

The unentangled triplet \T Ul z ) evolves similarly to 
\Tu T ,z}- Due to the continuous spin rotation induced by 
the s-o coupling, intermediate degrees of bunching and 
antibunching are possible as the modulation angle 8 so is 
varied. Note that a variety of entangled and unentan- 
gled electron pairs (along distinct quantization axes as 
well) can be generated according to the above states by 
sending electron pairs through a s-o active region. See 
also Ref. |6^ for entanglement generation using a Mach- 
Zchndcr interferometer. 
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where S^!^ a denotes the noise expectation value, Eq. 
(|46|l . for a pair with a spin-down electron in channel b 
of lead 1 and a spin-up electron in channel a of lead 2, 
I l*i)ba,z- Similarly, S^l^ a corresponds to the two oppo- 
site spin electrons in the channel a of the respective lead. 
The noise contribution of the pair state | []) a a,z can be 
determined from Eq. by setting 9 so = n which makes 
\T UTz ) L = -I H) aa , z , see C3, 



3. Coupled-channel case 

Here we assume that electron pairs are injected into 
channel a of leads 1 and 2 with energies near the crossing 
of the bands at k c , Figs. 1(b) and 2. The case with 
injection into both channels a and 6 is discussed at the 
end of this section. General expressions similar to Eqs. 
(|52|l and (|54|) can, in principle, be derived for this case. 
Here, however, we generalize the approach outlined at 
the end of the previous section to the the case of two 
s-o-coupled channels, which is both simpler and more 
intuitive. An additional heuristic derivation of the noise 
properties using simple number operators is given in the 
Appendix IC II 

a. Quantization axis along z. The idea is essentially 
the same as before: for instance, an injected unentangled 
spin-up triplet along the z direction, with one electron in 
channel a of lead 1 an the other in channel a of lead 2 



S. 



33:aa 



2f 
hp 



T(l-T). 



Note that the pair state | [^) aa ,z is distinguishable (spin 
down in lead 1 and spin up in lead 2, both in channel a), 
and therefore yields the classical shot noise. Similarly, for 
the distinguishable state | ||)fca,z we should have S. 



IT,* 
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S sl-L- Hence we find 
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^33 



(0 so ,0 d ) = ^-T(l-T)[l + sm 2 (d d /2)/2S £u 



cos (0 d /2)cos (9 so )6 £u 



(89) 



and, similarly, S^ 1 '* (6 so , d ) = S%?>' {6 s0 ,6 d ). The 
noise for the entangled triplet and the singlet are, 
respectively, 



$33" (fi'O, 0d) 
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l--(cos 2 (0d/2) + l)5 eil e a 
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and 



2e 



S$3(0 SO , 9 d ) = —T{l-T) [1 + cos {6 d /2) cos (9 so ) S eu£2 ] . 

(91) 

b. Quantization axis along y. We can again straight- 
forwardly find all the shot noise expressions for the quan- 
tization axis along y by following the same procedure as 
above. However, here a spin up (along the y direction) 
electron injected into the channel a of lead 1 evolves into 

|T)^ - ^{[(^ncos(^/2) e -^/ 2 



+ (l-iOe" W2 ]| T)a„ 
+ [(£ + i)cos(0 rf /2) 



+ e»»' 2 ]\l) a ,. 



+ (e-i)sm(6 d /2)e- ie °°/ 2 \ T) M 
+ ^-l)sin(V2)e-^ /2 |4) M }, 



(92) 



upon traversing the s-o region in lead 1, while a spin 
down electron (along the y) evolves according to 



■vl a,y 



= ^{[(l-C)cos(0 d /2)e-*W2 

+ [(£ - i) cos (0 d /2) e- l9 °°/ 2 - (i + e^° /2 ] | |) 0) , 

- (i + Osin(0 d /2)e- i( W 2 | T)m 

+ (^+l)sin(^/2)e-^°/ 2 | |) M }. (93) 



due to the weak coupling between channels a and b at 
the crossing. 

From the states ()94|l and (|95[) we can construct the 
electron pairs (entangled or not) that we are interested 
in, e.g., the singlet and the triplet states along the y 
direction. Let us consider a pair injected into only the 
channel a of leads 1 and 2. For the spin- up triplet we 
find 

\T Ur J L =e X p(i0 fl /2)|TT)aa,„, (96) 
while for the spin down 

\T Ul , y ) L = expH0 R /2)[cos(0 d /2)| U)aa, V 

- isin(0 d /2)| U>6a,»]. (97) 

The entangled triplet and singlet are found to be 

\S/T ey ) L = -^[exp(ife/2)| U)aa,y 

T exp(-^ H /2)cos(^/2)| |t)6a,j/] 

± -Liexp(-^ il /2)| n)ba, y , (98) 

where the upper (lower) sign corresponds to the singlet 
(triplet) state. The noise expression corresponding to the 
above states are 

(0 d ) = ^T(l — T) (1 — 5 ei , £2 ) , (99) 



For simplicity, we consider below the case with only 
the Rashba coupling, i.e., a/0 and (3 = 0. In this case, 
the above equations simplify to 



I /a,y 



a,y, 



(94) 



and 



| i) a ^ y L = exp(-i6 R /2) [cos(0 d /2) | l) a . y -iwn(0 d /2)\ ]) b J 

(95) 

where 6 so ^ 9r = 2maL/h 2 and | T)a,y and | |) a ,i/ are 
the eigenspinors of o~ y in channel a (similar definitions 
hold for channel b). Note that the spin- up and spin-down 
states in channel a of lead 1 evolve quite differently. In 
particular, the spin-up state remains in channel a and 
only picks up a overal phase, while the spin-down one 
acquires a spin-up component in channel b. This distinct 
evolution of | 1)a, y and | Day ca n be understood if we re- 
call that the incoming electron state | 1)a,y is essentially 
an eigenstate of the system at x = (far away from any 
level crossing) and hence evolves as such, as the electron 
traverses the s-o region. In contrast, the state | Da,y 
results from the incoming spin-down state in channel a 
being injected at the level crossing corresponding to the 
states | Da,y and | Db,y and thus evolves as a coherent 
superposition of these two states along lead 1. Hence, 
| Da y ls an eigenstate of our Hamiltonian in the approx- 
imation we considered in Eq. Q26[l. while | Day is n °tj 



S. Uv 
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and 



S 33 " y ^R^d) 
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T(l-T) [1 ± cos(0 d /2) cos(e R )5 £u 



(101) 

, Interestingly, the above results show that by measuring 
the noise in lead 3 along the y quantization axis one can 
distinguish all the electron pairs as they display distinct 
noise for non-zero s-o induced intcrchanncl mixing angles 
{9 d 7^ 0). Figure 3 shows the reduced Fano factor / = 
S33/ [2e 2 T(l - T)/hv\ for the singlet and triplets along 
the y axis as a function of 9 so and 9 d . It clearly shows 
that the singlet and triplets display distinct shot noise in 
a wide range of angles. Note that this result holds only 
for incoming electron pairs injected initially into channel 
a and with energies near the crossing. 



C. Effect of backscattering in the Rashba lead 

Here we examine how the presence of backscattering 
in the lead containing the spin-orbit interaction (lead 1) 
affects our results for the current fluctuations. A small 
backscattering can indeed be produced by the small band 
offset between the regions of the lead with and without 
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the Rashba interaction; see Appendix iBl where we calcu- 
late explicitly the transmission and reflection coefficients 
for a model quantum wire with spin-orbit interaction of 
both the Dresselhaus and Rahsba types. The situation 
is depicted in Fig. Ufa). 

We emphasize that here the "left" leads (1,4) and 
"right" leads (2,3) are not equivalent anymore because 
of the presence of the spin-orbit region in lead 1. This, 
as we shall see below, yields that S33 7^ S44. The case 
in which backscattering is present in both leads 1 and 2 
was first addressed in Ref. [63, and is discussed in the 
Appendix IC II within a simple heuristic picture [see Fig. 
H{b)]. For simplicity, we consider here only the case of 
uncoupled bands (Sec. II-C). 

We take into account the backscattering by adding a 
tunnel barrier with reflection amplitude a to lead 1, de- 
scribed by the scattering matrix 



Sback — 



(102) 



In order to combine Sback with the beam-splitter scatter- 
ing matrix s [Eq. (|47fl ] , we match the coefficients of the 
wavefunction going out of Sback with the corresponding 
ones going in s, and solve the resulting system of equa- 
tions eliminating the intermediate coefficients^. We ob- 
tain the modified scattering matrix 



a br bt 

t r 

br t ar 2 art 

bt r art at 2 



(103) 



The probability A = \a\ 2 of backscattering in lead 1 



renormalizes the transmission probabilities R 



and 



T = \t\ 2 by the factor \b\ 2 = 1 - A. It is interesting 
to note that the constraints of unitarity introduce some 
backscattering in the leads 3 and 4. These correspond to 
the possibility for particles injected into leads 3 or 4 to 
scatter back into leads 3 and 4 after being reflected by the 
tunnel barrier in lead 1. We note that the probabilities 
R and T are the ones of the original beam-splitter, and 
as such satisfy the normalization condition R = 1 — T. 

We proceed as before and introduce the spin rotation 
due to the Rashba coupling by multiplying S13 and S14 
by U so , Eqs. (|4*5|) and |39jl. The spin rotation has no 
effect on su even if the backscattering occurs after the 
Rashba region, because the spin of a particle with re- 
versed momentum rotates in the reversed direction. We 
can now use the general formulas Eqs. I|52|l and 1541) to 
evaluate the autocorrelation noise in lead 3. Introducing 
the renormalized probability R' = R(l — A), we find 



^33 ( 



,0) = j^[TR + R'(l- R') + 2TR' cos(0 so )tf £l , 



J 33 



hv 



(104) 



[TR + R'(l - R') - 2TR'5 eue2 ] , (105) 



S 3 P"(0 so ) = — [TR + R'(l - R') - 2TR' cos 2 (6 so /2)S eue2 ] 

(106) 

with the average current in lead 3 



h = V-(l-AR). 



(107) 



The factor 2TR found in the previous case without 
backscattering, Eqs. (|61J) ~ H63(1 . splits into two contribu- 
tions, TR + R'(l — R') with the renormalized probabil- 
ity R' . The backscattering adds a contribution related 
to the partition noise created by the tunneling barrier. 
Partition noise (shot noise) corresponds to fluctuations 
arising from the fact the barrier splits randomly the in- 
cident electron flow into transmitted and backscattered 
flows. On the other hand, the backscattering reduces the 
transmission probability R' for electrons injected in lead 
1, and therefore can decrease the "beam-splitter noise" 
proportional to R(l — R). For instance, the noise for 
entangled triplets with same energies is increased from 
zero (in the absence of backscattering) to a finite value. 
Similarly, the noise for electrons with different energies 
is increased in the range < A < max{2 — 1/T, 0}, with 
a maximal increase of (T — 1/2) 2 (e 2 /hv). However, the 
maximal value for the noise (obtained in the case of sin- 
glets with ei = &2 and 6 so = 2mr, n G Z) is not changed, 
as it corresponds to the maximal value reached in the 
case of bunching of bosonic-like particles. In the case 
of spin pairs injected along the Oy axis, the result for 
the singlet remains the same, while S 33 " y = S 33 z and 
S^(8) = Sl 3 (7T-9 so ). 

Figure [S] shows the auto-correlation noise 633 in the 
presence of backscattering in lead 1 , of the singlet S and 
for the unentanglcd triplet T UtZ in the case of equal en- 
ergies, ei = 62- We recall that the noise of the entangled 
triplet T e>z is independent of 9 so , as are S and T U)Z for 
different energies. We see that the backscattering can 
either increase or decrease the noise, depending on the 
rotation angle 6 so as well as the beam-splitter properties 
(defined by R). However, in all cases the backscattering 
reduces the visibility (amplitude) of the 9 so oscillations. 
In the absence of backscattering (A = 0), the bunching 
of singlets for 6 so — is maximal (S33 = e 2 /hv) for a 
symmetric beam-splitter (i.e., T = R = 0.5); this comes 
from the fact that the randomness is maximal in this 
configuration. In the presence of backscattering, the ef- 
fective transmission is decreased via R — > R'(l — A), so 
that one moves away from this symmetric point, the ran- 
domness of the beam-splitter is reduced, and the noise is 
suppressed. Similarly, the perfect anti-bunching seen at 
8 so = 7r (and at 6 so = for the triplets) no longer hap- 
pens because of the residual partition noise related to the 
backscattering and the noise is given by S§ 3 (8 so = n) = 

(8 so = 0) = S^t (ei = e 2 ) = (e 2 /hv)AR(l - AR) 
(see AppendixEl where we present a heuristic derivation 
of these results via simple calculations of expectation val- 
ues of number operators). 
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In the case of an asymmetric beam-splitter with a large 
value of R > 3/(4 - A) (e.g., R = 0.9), the noise for 
the singlets is increased by the backscattering for all 6 so . 
For 8 so = 0, one can understand this by noting that the 
backscattering reduces the value of R' which therefore 
gets closer to the symmetric value of 0.5 related to the 
maximal noise. On the other, in the case of small R (e.g., 
R = 0.1), the reduction of R' moves even further away 
from the maximal 0.5 value, and the noise is reduced 
for small 9 so . Around 6 so = tt it increases from zero, 
because of the additional partition noise. For the triplet, 
the maximal value is reached at 8 = tt, where electrons 
reach the beam-splitter with opposite spins, and hence 
behave classically. The value £33 = 0.5e 2 /hv at 6 so = tt 
for A = is modified in the presence of backscattering; 
it is then increased when R > 1/(2 — A). 

For the cross-correlations of the noise between leads 3 
and 4, we find 
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h v 



TR + T'R' + 2TR' cos(0)<5 £l , 
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TR + T'R' - 2TR'S C 



(108) 



(109) 



T„, 
°34 



(9) =-j- [TR + T'R' -2TR' cos 2 (6 /2)6 tl , t2 ] 



(110) 

with T' = T(l — A). These expressions are very similar 
to the ones for the auto-correlations, Eqs. (|108|) - (|110fl . 
The main difference is the negative sign of these correla- 
tions, which comes from the fact that S33 + S34 = when 
A = 0. The other difference is that the term R'(l - R') 
is replaced by R'T' ', which is easily understandable since 
the cross-correlation involves (713) (714) = R'T' for a parti- 
cle injected in lead 1, while the auto-correlation involves 
(n 3 )(l - (n 3 » = R'(l - R'), see Appendix© 



D. Spin-polarized case 

We now turn to the case of injection from a spin- 
polarized Fermi liquid lead with a continuous energy 
spectrum. We start from the standard expression for 
the noised, which corresponds to the continuous version 
of Eqs. (|44I46|) 

Safi = ^Rc JdeJ2 A^(a,e,e)A^(P,e,e) x 



■ySuv 



fju(e) [1 - fsv(e)} , 



(111) 



where f au (e) = {(J au a au ) and u = (a,n). Here, again, 
a =t,l denotes the spin components and n = a,b the 
band index. At zero temperature, for the scattering ma- 
trix of the beamsplitter with a two-band spin-orbit active 
region in one of its incoming arms (lead 1), Eqs. (|48|l and 



(14911 . we obtain for the shot noise in the outgoing leads 
(say, lead 3), 

- j JdeT{i-T) y, w^mWh^m- 

a a' rim 

ho> m (s)\ + han{e) [1 - f la > m (e)} }, (112) 

where U^ c n CT , m denote the matrix elements of \5f Q given 
in Eq. (g2J). 

1. Injection into one subband 

We consider a (non-equilibrium) spin-polarized injec- 
tion into leads 1 and 2 in subband a only, which we model 
with 
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(113) 
(114) 
(115) 



The degree of spin-polarized is controlled by the param- 
eter < p < 1. For full polarization p = 1 there is no 
voltage drop for spin-down electrons in Eq. (|114fl , which 
therefore do not contribute to transport. Note that here 
we only inject electrons in channel a, see Eqs. <|115H an d 
(|ll(j[l . We then obtain for the current in the outgoing 
leads 



h = h = ^V- 



1+P 



(117) 



which gets halved for full polarization, p = 1, as com- 
pared to the unpolarized case p = 0. For the auto- 
correlation noise we find 



2e 2 

S33 = — T(l-T) 
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\UTa;la\ 2 



P 
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(118) 



Note that S33 = S44 = Wc define the reduced 

Fano factor f p ee F/T(l - T) = S 33 /2eI 3 T(l - T). and 
using Eq. (|4^|l and |^| 2 = 1 we find 

, P ( 2 ^ 1 1 9 h \, 1 . 2 0d 
Jp = ^ I cos — hi 2 cos — cos e so 1 + - sin — . 

_ (H9) 
This agrees with Ref. |63 for (3 = and p = 1. Notice, 
however, that if p < 1 there is an additional correction, 
i.e., the last term in Eq. H119(l . which is independent of 
p and thus there is a small noise power even for the un- 
polarized case p — 0. This is related to the presence of 
the additional empty channel &, which allows for flucta- 
tions in the outgoing leads in the case of finite s-o in- 
duced interband coupling, 9d 7^ 0. In the strictly ID case 
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(i.e., single channel), this effect disappears as 6d = and 
Ecr. (|119|) reduces to the expression 

/ P = |(l-cos^ )-psin 2 ^. (120) 

This effect also vanishes in the case where both bands 
are occupied, studied in the next section. 



2. Injection into both leads 

We now consider the more realistic case of injection 
into both subbands a and b and model it with 
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where /ii j(T . a = Mi,a-,f> = ^i.a- We obtain the current 



2e 2 2 
h = h = -r-V—- , 
h 1+p 



(121) 
(122) 
(123) 

(124) 



which is twice the current in the case of injection into a 
single channel given by Ea. l|117fl . The Fano factor reads 



F = T(1-T)£ \Wn, 



(125) 



7i.rn—a,b 



Using Eq. (|4*2*|) . we obtain the reduced Fano factor f p 
F/T(l - T), 



f p = - I 1 - COS — COS 6 so 



(126) 



This vanishes in the case of unpolarized injection, p = 0, 
showing that the interband coupling in lead 1 in itself 
does not give rise to additional noise for two filled chan- 
nels. For an uncoupled channel (9 d — 0), f p yields the 
usual form valid for single channels, i.e. Eq. (|12L)fl . The 
above results were derived for a quantization axis along 
the z direction. For other directions, one can derive simi- 
lar but distinct formulas. Finally, we note that measuring 
such Fano factors enables one to quantify the degree of 
spin-polarization p of the reservoirs along different quan- 
tization axis. 



E. Coherent injection into multiple discrete states 

Here we generalize our previous results by studying 
the injection of electrons into leads with a discrete spec- 
trum. We show that even in the case of injection into 
many different discrete states of the lead, we can observe 
two-particle coherence, e.g., bunching and antibunching, 
and thus detect entangled states. We will identify six 
asymptotic regimes separated by the relative magnitude 



of the level spacing 8 in the leads, the energy mismatch 
A of the injected electrons, and the energy broadening 7 
of the injected electrons. It turns out that in four of the 
six regimes, we obtain (asymptotically) full two-particle 
interference. Only in two regimes, namely when the en- 
ergy mismatch A exceeds both S and 7, we obtain no 
two-particle interference. 

Electron injection into leads 1 and 2 is assumed to be 
coherent but with a finite width 7 in energy, such that 
several energy levels of the lead can be filled. This rep- 
resents a generalization of both the injection into single 
discrete levels as discussed in this paper and in earlier 
works&S^ and of recent investigations of the problem in 
the continuum limiljSi which will be a special case of the 
following discussion. 

We assume here injection into leads with equidistant 
levels e n = n5 + qd, where n — 0, ±1, ±2, . . . and < 
q < 1 is a fixed fractional offset. The assumption of 
equidistance simplifies our calculations because it allows 
us to perform the discrete sums explicitly (Appendix |EJ) . 
Also, we choose our discrete index n to run from —00 to 
+00, this is done for convenience and merely means that 
the spectrum is well described within an energy band 
a few times 7 wide. We do not expect our results and 
conclusions to be significantly altered in cases of non- 
equidistant energy levels in the leads. 

Injection of an electron with spin a =|,| into lead 
a and centered about the energy e is described by the 
creation operator 



E 

n— — 00 



g(e,e n )ai a (e n ) 7 



(127) 



where a^ aa (e n ) creates an electron with the sharp energy 
£„, cf. Eq. (|44|l and below. The weight function g will be 
assumed to have the Breit-Wigner form 



£ — £' 



I'J 



(128) 



with the normalization condition X)^L-oo lff( £ ' £ ')| 2 = 1- 
Unlike for the case of weak tunneling from, e.g., a quan- 
tum dot£, we assume here that the single-electron states 
arc filled with probability 1, but with uncertain energy 
This can be achieved with time-dependent tunnel bar- 
riers. Note that for widely spaced levels 6 ^> 7 and 
away from symmetric injection into two adjacent levels 
(q = 1/2), injection thus takes place into the nearest level 
with probability ~ 1. The two-particle injected states we 
arc interested in arc 



\S) \_ 1 



— (4 T (ei)4 i (e2)T4 i (Ei)4 T (e 2 )) , 



V2 

(129) 

the singlet and entangled triplet states with single- 
particle energies centered around e\ and £2 and smeared 
over a width 7. Using l|127|) . these states can be expressed 
in terms of the states with sharp energies Eq. (|5U)l as 

\S,T ei ) =^ e , E ,g(e 1 ,e> 1 )g(£ 2 ,e' 2 )\S,T ei ) e , E ,. Using the 
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normalization condition, we find that the average current 
in the outgoing leads of the beamsplitter is unaffected by 
the spread in energy, i.e., 1$ = I4 = —ejhv. For the Fano 
factor F = S^/2eI^, we find 



F S ,t. = T(1-T)(l±\h(e 1 ,s 2 )\ 2 
The discrete overlap function h is given by 



h{e ll e 2 ) 



g{e 2 ,e)*g{s 1 ,e). 



(130) 



from Eq. (|132|l . The limit 7 « k< A describes the situ- 
ation where the two electrons are injected into two differ- 
ent discrete states and thus \h\ 2 — > 0. In the other limit, 
7 « A « the two electrons are injected into the same 
discrete level and \h\ 2 = 1 - n 2 A 2 /3S 2 + 0(A 4 /S 4 , 1 2 ). 
The limit 7 <C S, A discussed here describes injection into 
single discrete levels discussed throughout this paper and 
in earlier works&SSiSi where £12 are multiples of i5 and 



(131) \hf 



is given by a Kronecker delta. 



The discrete summations required to evaluate the func- 
tion h are carried out in Appendix [E| and lead to the 
main result of this section (A = £2 — £1), 



N 2 = 



7 



cosh 2 (27T7/(5) - cos 2 (7rA/5) 



(A/2) 5 



sinh (2wy/5) 



(132) 



Note that \h\ 2 appears in Eq. I|130|) for the Fano factor 
and generalizes the Kronecker delta for the case of injec- 
tion into a single level. The function < \h(A, 5, j)\ 2 < 1 
determines the amount of two-particle interference that 
can be observed (see Fig. 01. Note that as a consequence 
of the normalization of g(e, £„), the overlap function \h\ 2 
is independent of the energy offset q (see Appendix lE|l . 
For \h\ 2 = 1, full bunching and antibunching of singlets 
and entangled triplets can be expected which is ideal for 
the purpose of discriminating these entangled states from 
uncntangled two-particle states. If \h\ 2 — > 0, then no in- 
terference, i.e., no bunching or antibunching will be ob- 
served because the two wavepackets centered around £1 
and £2 have no overlap in energy space. 

Since |/i(A, S, 7)| 2 depends on three parameters with 
the dimension of an energy, there are 3! = 6 different 
asymptotic regimes. We will now discuss all six cases, 
grouped into three sections with A, 7, and 6 as the small- 
est energy, respectively. We also note that since \h\ 2 is 
a dimensionlcss number, it is in fact only a function of 
two dimensionlcss energy ratios, e.g., A/S and 5/j, (see 
Fig. 0) as can easily be seen from Eq. l|132fl . 



1. Matching energies £1 = £2 (A <C *y,5) 

If A = £1 — £2 = 0, then \h\ 2 = 1 for arbitrary values of 
the other three parameters S, 7, and q. In other words, if 
the two electrons are injected with energy distributions 
whose centers coincide, then the interference is always 
maximal. This result persists for finite 7 and 6 as long 
as A < 7,(5; in this case, the correction is <3(A 2 ). For 
A < 7 < S, we find \h\ 2 ~ 1 - n 2 A 2 /35 2 . In the case 
A < 5 < 7, the correction is \h\ 2 ~ 1 - A 2 /A^ 2 . 



2. Sharp energies 7 <C S, A 
In the limit of sharp energies, 7 -C 5, A, we obtain 

N 2 = 4x? sin 2 t4 + 0( 7 2 ) (133) 



3. Continuum limit S <C A, 7 

In the continuum limit, i.e., the asymptotic case where 
the single-electron level spacing in the leads becomes 
small compared to all other relevant energy scales, we 
find a LorcntzianSi 



\h\ s 



T 



(A/2) 2 +7 2 ' 



(134) 



with the usual special limits \h\ 2 
and \h\ 2 — > 1 for 6 < A <C 7. 



for 8 < 7 < A 



4- Discussion in terms of largest energy scale 

Above, we have discussed all six asymptotic cases, 
grouped according to the smallest energy scale. Alter- 
natively, we can discuss the asymptotic regimes charac- 
terized by the largest energy scale. In the detuned case 
A 3> <5, 7, we never see two-particle interference since 
\h\ 2 ~ (4 7 2 /A 2 )/(A;<5,7) -> where /(A; S, 7) is a func- 
tion which is bounded for fixed S and 7. In the case of a 
wide distribution, as e.g. effected by a fast injection into 
the lead with injection time ft/7, we are in the regime 
7 ^> (5, A, and we always find that \h\ 2 exponentially ap- 
proaches 1. In the single-level case S S> 7, A, there is 
only one level to fill and we obtain \h\ 2 — > 1 irrespective 
of the relative magnitude of A and 7. 



IV. SUMMARY 

We have carried out a thorough study of current and 
noise for spin-polarized and spin-entangled electrons in 
a beam-splitter geometry, including a local spin-orbit in- 
teraction (Rashba and Dresselhaus) in one of the incom- 
ing arms. We have considered incoming leads with one 
or two channels, as well as backscattering effects. The 
channels can be coupled via the s-o interaction for in- 
coming energies near the band-crossing. We have found 
that the spin-orbit interaction is a useful mechanism to 
coherently rotate spin states. Such rotation can be used 
to modulate noise signals, thus providing unique signa- 
tures of spin-polarization and spin-entanglement. 

For spin-polarized electrons, noise measurements can 
give a direct measure of the degree of polarization along 
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different directions. For electron pairs, the coupling be- 
tween the channels can play an important role. For pairs 
with incoming energies near the band-crossing injected 
into one of the channels, we find an additional modula- 
tion due to the coherent transfer of electrons between the 
two channels. In this case, noise measurement allows us 
to distinguish all the different triplets states defined along 
the y direction, in addition to the singlet. Furthermore, 
for equal strengths the combined effect of the Rashba and 
Dresselhaus interactions can partially cancel out. In this 
case, the spin and orbital degrees of freedom are sepa- 
rable, the intcrband coupling essentially disappears, and 
the propagation of spin states is robust against scattering 
off non-magnetic impurities. 

We have also considered the influence of backscattering 
in the beam-splitter with a single channel. The main ef- 
fect is an additional contribution related to the partition 
noise due to the tunnel barrier describing the backscat- 
tering. This reduces the visibility of the oscillations in 
the shot-noise as the s-o rotation angle 9 so is varied. It 
also reduces the maximal noise value found for perfect 
antibunching of singlets. 

We have generalized earlier results for the shot noise 
of entangled electrons by allowing the injection of wave 
packets, i.e., coherent superpositions of discrete momen- 
tum eigenstates (plane waves). We have found a general 
analytical formula for the two-particle interference visi- 
bility \h\ 2 in terms of all three relevant energy scales A, 
7, and S. Our new result contains and generalizes both 
the discrete single-level case and the continuum case. 

Finally, we have developed a simple heuristic picture 
for the noise based on number operators in the different 
leads and the relevant transmission and reflection proba- 
bility amplitudes in the beamsplitter. Within this picture 
we can more intuitively re-derive some of the formulas for 
the noise - previously derived within the rigorous scat- 
tering formalism (Sec. Ill) - in the presence of spin-orbit 
interaction and backscattering in the incoming leads. 



APPENDIX A: BOUNDARY CONDITIONS AT 

x — L 



Here we show in some detail that both the wave func- 
tion and the velocity operator acting on it are continuous 
at the exit (x = L) of the s-o region in lead 1. 

1. Continuity of ty(x,y). 

The continuity of the wavefunction at x = L is trivially 
satisfied by the two-channel state 



cos (9 d /2) e-*W2 + e *„/2 \ l i{ko+kso) 



£[cos(0 d /2)e- i0 «<»/2 _ e *W2] j 2 
^sin(0 d /2)e-^°/ 2 ) 2 6 Myh 



describing the electron state within the s-o region, < 
x < L, and the state 



(Al) 



(A2) 



valid for x > L in lead 1, if we choose A, B, C, and D, 
equal to the corresponding components of ty(L,y). 



2. Continuity of the current flow. 

The continuity of the (non-diagonal) velocity 
operator— 
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acting on the wavefuctions at x = L 

v S o^(x,y)\ x ^L- = v so ®(x,y)\ x ^ L +, 



(A4) 



assures current conservation. The left-hand side of (|A4(I 
yields 



z} so *(x,t/)| a 
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£(fc 2 -fc so )e^ L 
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FIG. 3: Fano factors as a function of Or and Od for triplet 
(a) and singlet (b) pairs (defined along the y axis) with injec- 
tion energies near the channel anticrossings. We assume the 
injected pairs to be initially in channel a of the non-Rashba 
region in lead 1. In addition to the usual Rashba-induced 
spin rotation Or, the spin-orbit interaction induces a further 
modulation on the Fano factors via the coherent transfer of 
electrons between the channels (mixing angle 9d)- Interest- 
ingly, the triplets display distinct Rashba s-o modulations: 
the unentangled spin-up pair T U]V is not sensitive to spin- 
orbit effects and shows full anti-bunching, the unentangled 
spin-down triplet T UlV depends on only 9d and oscillates be- 
tween the full anti-bunching and the classical value for dis- 
tinguishable pairs, the entangled triplet T By displays sizable 
oscillations between full bunching and anti-bunching as Od and 
Or are varied. For the sake of clarity, we have omitted in (a) 
the Or dependence of Fano factors corresponding to T u ^ y and 
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FIG. 4: Two cases for a beam-splitter with backscattering. 
(a) Back-scattering in arm 1 only, induced by the s-o region, 
with a the backscattering amplitude, and r and t the original 
amplitudes of the beam-splitter, (b) Backscattering present 
in all leads (amplitude a), with also cross-backscattering (c). 
Here r and i are the modified amplitudes satisfying the nor- 
malization | a | 2 + |c| 2 + \r\ 2 + \t\ 2 = 1. This case will be 
discussed in the Appendix IC II 
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where we used k2 — k so = k c + fc so . On the other hand, 
it is straightforward to show that 
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2m 
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Again, assuming A « fc c ~ we can drop the terms 
proportional to A/ (k c ± fc so ) in Eq. l|A6p thus arriving at 
the desired equality v so ^f(x,y)\ x ^ L - = v so $(x 1 y)\ x ^, L + 
which assures current conservation. 



APPENDIX B: TRANSPORT THROUGH 
QUANTUM WIRE AT \a\ = /3 

As an example of a backscattering mechanism in quasi- 
1D channels, we calculate here the transmission and re- 
flection coefficients for a quantum wire of length a in the 
presence of spin-orbit coupling of both the Rashba and 
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Rashba angle 9 Rashba angle 9 




FIG. 7: The function \h(A,8,-y)\ 2 as expressed in Eq. IT321 . 
plotted versus the dimensionless quantities A/S and (5/7, 
where A = £2 — £1 denotes the mean energy difference be- 
tween the two injected electrons, 7 the width of their energy 
distributions, and 8 the level spacing in the leads. The case 
of matching energies E\ ~ £2 or A <C 8, 7 (Sec. 1111 E 11 cor- 
responds to the edge indicated as "small A"; here we find 
\h\ 2 = 1, irrespective of the ratio 7/5. The case of a broad 
energy distribution (fast injection) 7 ^> A, S (Sec. 1111 E 41 cor- 
responds to the edge indicated as "large 7"; along this edge 
\h\ 2 = 1, regardless of A/S. The red arrow follows a line 
of constant 7 and A with variable 8. For large 8, we find 
again \h\ 2 = 1 for all A/7, whereas for small 8 (continuum 
limit, Sec. llIIE3L the limit \h\ 2 = or \h\ 2 = 1 is reached, 
depending on whether A/7 > 1 or A/7 <C 1. 



FIG. 5: Noise in the presence of backscattering in the lead 
with s-o interaction, for backscattering probabilities given by 
A = 0, 0.1, and 0.4. The transmission probabilities are R = 
1 - T = 0.5 (a,b); R = 0.9 (c,d); and R = 0.1 (e,f). The left 
(right) column shows the case of the singlet S (unentangled 
triplet T u ,z) with ei = £2- 





FIG. 6: Injection of electrons into equidistant discrete lead 
states n = 0, ±1, ±2, . . . with spacing 8 at center energies £1,2 
separated by A = £2 — £1 with distributions <?(£i,2, £) of width 
7. The distributions g(ei : 2,£) are not drawn to scale, as their 
normalizations depend on q. 



Drcssclhaus types with equal strengths |a| = (3. The wire 
is attached to leads which are modeled as semi-infinite 
quantum wires without spin-orbit interaction. Again we 
choose the axis of the wire along the x-direction. The 
central region < x < a is characterized by an effec- 
tive mass TO2 and a confining potential V^y), whereas 
the leads (x < and x > a) have effective mass mi and 
confining potential V\{y) 1 in general mi ^ 7712, V\ 7^ V2. 
For dcfinitcncss let us consider the case a = ~(3. Here 
the spin is independent of the momentum along the wire 
thus being a good quantum number. The problem of 
matching the wave functions at the interfaces separates 
then into these two spin directions. Without loss of gen- 
erality, we will concentrate on the +- state as defined in 
Eq. Moreover, to analyze ballistic transport through 
this arrangement we will use the approximation that the 
reflected and transmitted electrons have the same sub- 
band index as the incoming ones. Thus, the incoming 
and reflected parts of the wave function at x < read 



Ae 



— iki. 



(Bi) 



and for the transmitted part at x > a we have 

i>(x,y)=Ce ikx c!> 1 n (y). (B2) 
Here k is the wave vector of the incoming particle along 



the wire axis with energy e 



(Hk) /2mi, where 



e n is the subband energy of the transverse wave function 
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ifin(y) according to the potential V\{y). In the above 
equations, A is the amplitude of the incoming wave, and 
A and C are the amplitudes of tne reflected and trans- 
mitted wave, respectively. Note that the above form of 
the transmitted and reflected part of the wave function 
are a restricted ansatz containing the aforementioned ap- 
proximation that the subband index is the same as in the 
incoming part. For the wave function in the region with 
spin-orbit coupling (0 < x < a) we use the general ansatz 



with 

= — \f2am-2. ± 
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and e 2 is the subband energy of the transverse wave func- 
tion ip 2 in the absence of spin-orbit coupling. The above 
ansatz contains all subbnads, and the corresponding wave 
vectors q± will acquire imaginary parts for large enough 
subband energies. However, since this region of the setup 
is finite (0 < x < a), the wave functions remain normal- 
izable. 

The continuity conditions on the wave functions at the 
interface x = lead to 
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and at 



one has 
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Multiplying the above by (4> 2 (y))*e- tV2mn2y and inte- 



grating over the transverse direction y gives 
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k{X-A)S jn = nJ!'. 'hH' 



— kCe ika S jn = Tl 3 Bie lq + a - mBLe 
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where we have defined the overlap integrals 
S jn = J dy{<j> 2 (y)ye + ^^y^ n {y). 
Eliminating the quantities A and C yields 
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or, solving for B±, 
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Moreover, multipying Eqs. I|B5|I - (|B9|1 by {4>n(y))* and 
integrating over y gives 



X + A = E [( B + + B -) S Jn 
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Inserting the above expressions for B± into Eq. (|B18(1 
yields the following result for the reflection amplitude, 
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and from Eq. (|B20|) one finds for the transmission ampli- 
tude 
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It is worthwhile to note that using Eqs. (|B19|l . ljB21|l in- 
stead of Eqs. l|B18(l . (|B20p leads to an identical expression 
for the transmission amplitude and to an equivalent re- 
sult for the reflection amplitude, 
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This is indeed the same as Eq. (|B22|) as one can see us- 
ing J^j \Sjn\ 2 = l- Thus, Eqs. (1b1T1 - |1bT1 (from which 
(|B17|) was obtained) are consistent with Eqs. (|B18|) - 
(|B21|I . This is a nontrivial finding since our original 
ansatz for the wave function was a restricted one con- 
taining an approximation, and it is not a priori clear 
that such an ansatz would lead to a consistent system of 
equations. 

APPENDIX C: HEURISTIC PICTURE OF THE 
NOISE IN A BEAM-SPLITTER 

One can obtain a simple physical picture of the ex- 
pressions for the current noise by considering the fluc- 
tuations of the number operator in the outgoing leads 
when one injects a pair of particles in leads 1 and 2 (i.e., 
(ni) = (712) = 1). For instance, 

S 33 (j^j ~ (An§) = (n§) - (n 3 )\ 



^ U \Jw) ~ (A?i 3 An 4 ) = (n 3 n 4 ) - (n 3 )(n 4 ). 

By considering classical, Fermi or Bose particles we can 
derive formulas which are in direct correspondence with 
the results for the current fluctations for electron pairs 
with spin previously obtained in Sec. Ill using the rig- 
orous scattering formalism. For clarity we discuss sep- 
arately the different configurations, illustrated in Fig. 
We first consider a symmetric beam-splitter with 



backscattering and no local s-o effect in lead 1; see Fig. 
0Jb). We then consider the case of backscattering in the 
s-o lead 1 only, shown in Fig. 0fa). Finally, we return to 
the s-o rotation in a beam-splitter with no backscatter- 
ing. 



1. Backscattering at the beam-splitter 

In realistic experimental the beam-splitter is not per- 
fect, and can actually contain a significant amount of 
backscattering in all input leads (i.e., sn,S22 7^ 0), as 
well as "cross-backscattering" between the input leads 
(S12 =/= 0). The most symmetric scattering matrix corre- 
sponding to this situation is 

/ a c r t \ 

c a t r I 

r t a c I ' 
\ t r c a ) 

where all backscattering amplitudes are the same s\\ = 
S22 = S33 = S44 = a and so are the cross-backscattering 
amplitudes S12 = S34 = c. Defining A = |a| 2 ,C = 
|c| 2 ,i? = \r\ 2 , and T = \t\ 2 , the unitarity of s imposes 
A + C + R + T = 1 and c = —ar/t, so that a and c 
are not independent. To have independent a and c, one 
must drop the requirment of symmetry between input 
and ouput, as is considered in section ITlI CI and IC 21 for 
the case of backscattering in the s-o lead only. We now 
calculate the expectation values of the number opera- 
tors in lead 3 and 4, by simply considering the probabili- 
ties for the different scattering configurations. Let P(3), 
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P(3,3), and P(3,4) denote the probabilities of finding 
one and two electrons in lead 3, and one electron in each 
lead 3 and 4, respectively. We need to determine (n 3 ) = 
P(3) + 2P(3,3), (n§) = P(3) + 4P(3,3) = (n 3 )+2P(3,3), 
and (713724) = P(3, 4). 

For classical particles, we have (n 3 )c = | S13 1 2 ( | S21 1 2 + 

M 2 + | S24 | 2 ) + |S 23 | 2 (|S11| 2 + M 2 + |S14| 2 ) + 

2|si 3 | 2 |s 23 | 2 = R + T, (n|) c = R + T + 2RT, and 
(ti 3 tm)c = |si3| 2 |s 24 | 2 + |su| 2 |s 23 | 2 = R 2 +T 2 . We find 
the auto- and cross-correlations 

(An 2 ) c = T(l-T)+R(l-R) 

= (A + C)(R + T) + 2RT, (CI) 



(An 3 An 4 )c = -2RT. (C2) 

As we shall see in App. IC 21 below, these results cor- 
respond to the current noise S 33 and S34 for electrons 
behaving classically (i.e., with different energies and/or 
with opposite spins). The result <|C1[) consists of the sum 
of two terms corresponding to the partition noise for elec- 
trons coming from lead 2 and 1, respectively. One can 
simply add these contributions because the classical par- 
ticles are independent. In the second equality we can 
recognize a partition noise term (A + C)(R + T) in addi- 
tion to the usual "beam-splitter noise" ~ 2RT. 

The situation is different for quantum particles obey- 
ing Fermi or Bose statistics, in which case one must first 
add or substract the amplitudes before building the prob- 
abilities for indistiguishable events. One can satisfy uni- 
tarity by choosing, for convenience, 3t(r*t) = 5R(a*c) = 0, 
which yields \r 2 ±t 2 \ = R + T and \a 2 ± c 2 | = A + C. 
For spinless fermions, one has P(3, 3) = and (n 3 ) F = 

|S1 3 S24 - Sl 4 S2 3 | 2 + IS13S21 - «U S 23 | 2 + | Si 3 S 22 ~ Sl2S 23 | 2 = 

R + T, (nl) F = (n 3 ) F = R + T, and (n 3 n 4 ) F = 
IS13S24 — si4S 23 | 2 = (R+ T) 2 . The correlations read 

(An 2 ) p = T(l — T) + R(l — R) — 2RT 
= (A + C)(R + T), 

(An 3 An 4 )^ = 0. 

For the autocorrelation (An 2 ) F , we see that the zero 
value found in the absence of backscattering becomes fi- 
nite, which is a consequence of the partition noise cre- 
ated by the additionnal backscattering channels. One 
can obtain this result by substracting from the classical 
result the forbidden case with two electrons in lead 3, 
(An 2 ) p = (A7i|)c — 2P(3, 3). The cross-correlations, on 
the other hand, remains zero. These results for fermionic 
particles correspond to electrons with equal energies in a 
triplet state. 

For spinless bosons, we must double P(3,3) thus ob- 
taining (n 3 ) B = R + T, (n§)s = R + T + ART, and 
(113114} b = (R — T) 2 . The correlations are 

(An 2 ) b = T(l — T) + P(l — R) + 2RT 

= {A + C){R + T)+ART, (C3) 



(An 3 An 4 ) B = -4PT, (C4) 

and correspond to electrons in a singlet pair with equal 
energies. We recognize again the sum of a partition-type 
noise and the beam-splitter noise 4PT. 

We can establish a connection between the above re- 
sults and those of Ref. in which the case with no 
cross-backscattering si 2 = was studied. In that work, 
the back-scattering was introduced in the same way as 
in Section IIII CI namely, by taking the probabilities R 
and T from the original beam-splitter and adding a tun- 
nel barrier with reflection probability Rb- Since the 
cross-backscattering does not play a direct role for the 
noise in the output leads, one can identify A + C — > Rb, 
R -> P(l - R B ), and T — ► T(l — R B ), which establishes 
the equivalence of Eq. lfCl|) . lfU4)l and Eqs. (5), (6) of Ref. 

E3 



2. Backscattering in one lead 

We consider here the case of backscattering in lead 1 
only, which was discussed in section ITlI CI The problem 
is no longer symmetric and the transmission probabilities 
T and R are not equivalent anymore. Sec Appendix B 
where we calculate explicitly the transmission and reflec- 
tion coefficients for a model quantum wire with spin-orbit 
interaction of both the Dresselhaus and Rahsba types. 
We recall that R and T here are the original quantities 
before adding the backscattering channel, and therefore 
satisfy T = 1 — R. We calculate the fluctuations of the 
number operator in lead 3 n3 by the same procedure as 
in Appendix IC II 

For classical particles we find (n 3 )c = 1 — AR, (n 2 )c = 
1-AR + 2TR, which yields 

(An 2 ) c = TR + R'(l - R). (C5) 

This result corresponds to the current noise expressions 
(|104|l . (|105|l . and (|1L)6[> with ei ^ e 2 . Hence, in this con- 
figuration electrons with distinct energies are not affected 
by the (anti-)symmetrization related to their spin state, 
and behave, effectively, like classical particles. 

For fermions, one finds (n 3 ) F — (n 2 ) F = 1 — AR, which 
gives 

(An§) F = TR + R (1 - R) - 2TR'. (C6) 

This result corresponds to the case of triplets l|l(J5|) , l|lUt)|) 
with ei = e 2 and 9 so = 0. 

For bosons we can proceed likewise to find (u 3 )b = 
1 - AR, (n 2 3 ) B = l- AR + ATR'. The result 

(An|) s =TR + R'(1-R') + 2TR\ (C7) 

corresponds to the case of singlets with ej = e 2 and 6 so = 
0. This can also be found by doubling the probability to 
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have electrons in the same lead 3, (An§).B = (An|)c 
2P(3,3). 

We can also rewrite the results above as 
(Anl) c = AR(1 - AR) + 2TR', 



(Anj) F = AR(1 - AR) = (n 3 )(l - (n 3 », 



(Arig) b = AR{1 - AR) + ATR' . 

This shows how the backscattcring induces a partition 
noise ~ (?i3)(l — ("3)), in addition to renormalizing the 
beam-splitter noise ~ TR' via R' = R(l — A). 



3. Spin-orbit rotation 



P(3) + 4:P(3,3) = {n 2 3 ) A +A 
l + 2TR(l + cos9). 



S13S23- 



V2 



(CIO) 



The factor (1 + e* e )/V2 interpolates between the Pauli 
exclusion principle and the bosonic bunching occuring for 
two particles in the same lead 3. 



APPENDIX D: EVALUATION OF THE NOISE 
MATRIX ELEMENT 



We can incorporate some of the effects of s-o induced 
rotation within the heuristic scheme presented above. 
We neglect backscattcring for simplicity. Considering 
the result for the singlet [Eq. ijpTj l] with ei = 
S§ 3 = (e 2 jhv) 2TR(1 + cos#), we see that the s-o angle 
8 so "interpolates" from the bosonic behavior (9 so = 0) 
to the fermionic behavior (Q so = ir). For the angle 
6 so = 7r/2, we actually recover the same result as in the 
classical case (ei 7^ €2) - although one must emphasize 
that any locally-rotated state is still maximally entan- 
gled, and therefore is not classical. In order to consider 
particles that effectively behave in an intermediate way 
between bosons and fermions, we consider particles that 
follow intermediate statistics, i.e., anyons that acquire a 
finite phase e z9 upon antisymmctrization (with 6 so = 
for bosons and 9 so = -k for fermions). Indeed, one can 
easily recover the formula 



(A? 



■3/A 



2TR{\ 



cos ( 



(C8) 



corresponding to a Rashba-rotated singlet pair with ei 
€2 [Eq. H6ip] by the following calculation for anyons 



(n 3 ) A = P(3) + 2P(3,3) = |si 3 S24 + e l %4S23| 



+2 



1 



S13S23- 



V2 



R + T=l, 



(C9) 



Here we sketch the derivation of general formulas for 
the relevant matrix elements in the noise calculation for 
injected electron pairs. From the noise definition [Eq. 
(|46|l ] it is clear that we need to evaluate objets of the 
form 



and 



D = (0|a ft<7(> (e M )o I/jeri/ (e 1/ )a^ o .(e) X 

x ap.a* (e')<4,<r 7 ( £ -t) a l^ ( s n)\°) ( D1 ) 



Q = (0|a Al , (Tfl (e /i )a 1/ , (T „(£ 1/ )a^ (7 (£)a 1 g, (T '(e') x 

x 4, a „ (e")a p/ ,^ (£">t iCT7 ( £7 )at ^ (e n ) |0p2) 



Wick's theorem tells us that we can express the above 
matrix elements in terms of all possible pairings of the 
fermionic operators. Here the possible non-zero pairings 
involve only one creation and one destruction operator. 
For the first case we have: 



D = +(0|at !CT (£)^ 1 ^(£')|0)(0|a,,^(£,)at i ^(£ 7 )|0)(0|a M ,^(£ /1 )at )CTj] (^)|0) 
-(0|at !CT (£)a^ lCT ,(£0|0)(0|a,, CT „(£,)at i(rr) (£,)|0)(0|a^ CTM (£^)a^(£ 7 )|0) 
+(O|a,, CT „(£,)at !CT (£)|O)(O|a0, CT ,(£Oa 7 ,^(£ 7 )|O)(O|a A1 , (r(1 (£ M )at !(7t) (£^)|O) 

-(0|a^(£,)ai iff (£)|0)(0|^, ff Ke')«U,,(^)l )( l^,^(^)«^ 7 (e7)|0) 
+(0|a^,^(£^)a^ CT (£)|0)(0|a, 1 ^(£,)a 7 ,^(£ 7 )|0)(0|a^ < (£ / )at )ffi; (£ n )|0) 
-(O|a M , CTfl (£ M )at^( £ )|o)(O| a0 ,^(£Oa 7 ,^(£ 7 )|O)(O| a ,,^(£,) a t^(^)|O). (D3) 
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The first two terms are zero since they both have a de- Similarly, we find 
struction operator acting on the vacuum state. Hence, 
we have 

-8p v 8 m 8 (T > a . ri 8 ( r^ (Trf 8 £ ' ev 8 eiMe ^) + 

8 t _ La 8a l ,a8e l ,e{8 vl 8j3 ri 5 c , va ^ 8 a > ar] 8 El/Ey 8 E ' Er} 

-8j3^8 1/ n8 cr '^8 (Tl/Crv 8 e ' ey 8 eu e n )- (D4) 
I 



Q = +5 a pS a ' ! 3'8 aa 'S a '' a >>>5 eE '8 E '' E >>>[8 lJ1 5 m 5 (JvCr _ Y 5 a ^ ari 8 evei 8 EfiEri 
~S a ip'S m S cr " a '"S crtirTv 5 E " E '"S EfiEv [Sv 1 Sp a S (7i/ay 5 cr ' rT S EiyE ^S E ' E 

+8 a 'j3'8 ll ^8 a ii (T "i <5cr M cr T S E » E '» 8 E ^ Ei [S url Sj3 a S auarl S a > a S EijE7i 5, 

—8 a 'i3'8 l j, a 8a" (T '"8 afia 8 e " e '"8 efie [8 V n80y8a u (T V 8 (T 'a y 8 £t/ev 8 l 

—8 va 8 m 8o v a8c ril a ri 8 £v£ 8 £fl£ri [8p 1 8p' a ' 8 <T 'a 1 8 a '"a" S E ' Ey 5, 
+ 8ua8 fll Sa v cr8a tl cr y 8 EijE S E ^ Ey [8pr,8p' a ' 8 a ' a v 8 a "' a" 8 B > 
—8va>8 l ict'8o u a8o tl ,tT"8e u s8 £liE " [Sp rl Sp' 1 5 a r (Jr] S ( j>/> cry 5 E ' Er] 5 i 
+8 l j, a 8 vv Sa fl a8a 1 ,a n 8 efie S Sl/en [8pr r S a t crj 8 £ ' £ ^8p' a '8cr"'a"8 l 

—8 iia 8v 1 8 atilT 8 lTu < 7l 8 £liS 8 E „ £ ^ [5p v 8 a ' C , ri 5 E ' Ev 8p' a '5 a '" Cr ii5, 

+8 l j, a S va >S aiir7 8 C r ua >>8 £fl£ S El/£ ri [8p v S a ' r7n 8 £ t £v Sp' 1 S a '"a y 8 l 



e's 
s's-y 
e"'e" 
e"'e" 
e f/f s y 
e'" e" 
e'" e" 

€ f ' f €-y 



8fj,y8 vr i8 crilCr ^ 8 av(7ll 8 EiiE _ i 8 EijEv \ 
8va8p 1 8< 7vC ,8 a ' a _ i 8 EvE 5 E i £y ] 

8va8pr)8cr v a8(T>a n 8 £l/£ 8 £ i £f) \ 

8vydpnScr v <j^ 8<j>a 7] 5 EuE _ i 8 £ > £ri ] 

8pa'8p' 1 S a ' cr "S a '" (Jy 8 E ' E "S E '" E _ i ] 
8 pa 1 8 a ' a" 8 E ' E " 8 pi ^8 a i'i 0r] 8 E »' Erl ] 
8p 1 8pi 71 8 a 'a 1 8cr>»v v 8 E ' El S E "' Erl ] 
8p a '8 a 'a" 8 E ' E " 8p' 7 8 a "' a ^ 8 £ "' £y ] 
8pa'8a'a"8 E ' E "8p' v 8 a >" a71 8 E "' Ev ] 
8p 1 8 a i a , l 8 E i E , f Sp'^Sa'i'a^ S E "' Erl ] . 



(D5) 



With the help of Eqs. (|D4(I and 1D5(I we can system- 
atically determine all the relevant matrix elements ap- 
pearing in the noise calculation for a particular type of 
injected electron pair (singlet, triplets, Bell states, etc.). 



particle spectrum of the leads, e = E\ + (g + n)8, where 
n is an integer, q a real number between and 1, and 8 
is the level spacing of the leads, see Fig. [5] We obtain 
h(q) = go(q — — )go(q)ho(q) with the discrete sums 



APPENDIX E: DISCRETE SUMS FOR h(s 1 ,e 2 ) 

We choose e\ as our reference energy and define A = 
£2 — E\. We further assume an equidistant discrete singlc- 



1 _ -2tt sin Tr(q + (cos ?r(g - 3i|) - cos ?r(g + i% - 2f )) 

~~ n =^oo ( nS + 1 S - A + h){nS + q8- i 7 ) ~ 8{A - 2i>y) (cosh27rJ - cos 2nq) (cosh 2ir} - cos27r(g - £)) ' 

(El) 

1 1 7T sinh27T^ 

9o(q) 2 _ ( n 8 + Q8) 2 + 7 2 7<5 cosh 2ttj — cos 2irq ' 



Taking the modulus squared of the complex function h, 
we Obtain |%)| 2 = g (q) 2 g (q- A/S) 2 \h (q)\ 2 = AB(q), 
with a manifestly g-independent factor 



.4 



(A/2) 2 + 7 2 2sinh 2 2tt 7 /(5' 



(E3) 



and an apparently g-dependent factor 

| cos7r(g — 3i-f/S) — cos7r(g — 2A/<5 + i"f/8)\ 2 



B(q) 



cosh27r 7 /<5 — cos27r(g — A/8) 



(E4) 



Inspection of Eq. (|E4(1 shows that the dependence of B(q) 
on the offset g drops out completely, and we arc left with 
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B = B{q) = cosh(47r 7 /(5)-cos(27rA/5) = cosh 2 (2n-//5) - 
cos 2 (it A/ S), where we have used cos 2 a; = (1 + cosx)/2 
and cosh 2 x = (1 + cosh a;) /2. Combining A and B, we 
finally obtain Eq. I|132|) . We note that the cancellation 



of q from the interference function \h\ 2 is a consequence 
of our normalization condition, i.e., that electrons are 
injected with unit probability. 
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